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Hurl r>'  the  paut  1/  y**arc  of  research  on  the  subject 
"Theoretical  and  experimental  studies  of  VLF  and  LF 
waves''  quite  a  number  of  scientific  conclusions  could 
be  drawn  that  are  reflected  in  more  than  hundred 
papers  and  were  read  at  several  scientific  confer¬ 
ences  and  meetings. 

Under  the  VLF  Research  Contracts  nine  comprehensive 
thesises  were  completed  up  to  now,  four  more  are 
in  preparation. 

experimental  measurements  and  investigations  were 
conducted  in  the  home  mines  of  Grosskogel/St .  Gertraudi, 
Falkenstein/Schwaz ,  Lafatsch/Scharnitz ,  and  salt  mine 
Hall  (all  of  them  in  the  Tyrol) ,  as  well  as  in  Blei- 
berg/Carinthia,  Mitterberg/Salzburq  and  in  some  German 
and  Dutch  mines. 

In  addition  to  commercial  apparatuses  and  measuring 
instruments  we  have  designed  and  built  a  great  number 
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efficiency  for  our  specific  purposes.  We  may  especially 
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St.  Gertraudi,  e.g.  the  fully  conditioned  mine  labo¬ 
ratory,  the  flood able  gallery  with  the  pumping  equip¬ 
ment  for  subsurface  experiments,  and  the  big  frame 
antenna,  that  could  be  constructed  despite  of  our 
limited  budgetary  funds  thanks  to  the  generous  assist¬ 
ance  of  our  home  industries. 

The  general  rise  in  prices,  the  shortening  of  budgetary 
funds  in  connection  with  the  financial  situation  of 
the  USA  and  last  not  least  the  devaluation  of  the 


dollar  by  more  than  20  %  In  the  course  of  the  last 
two  years  sometimes  brought  about  very  serious 
problems  that  could  only  be  solved  by  the  working 
morale  of  every  individual  member  of  the  VLF  team. 

We  are  indebted  to  the  director  of  the  Montanwerke 
Brixlegg,  Dipl.-Ing.  Paul  KETTNER,  owner  of  the  mine 
of  St.  Gertraudi,  for  being  allowed  to  conduct  our 
measurements  there. 
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Univ.-Prof .  Dr.  Rudolf  STEINMAURER,  Head  of  the 
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members  of  the  VLF-Project  for  the  good  and  success** 
ful  co-operation . 


Norbert  Nessler 


Ill 


VLF  Tean.  1972 


Contractor:  Dr.  Norbert  NESSLER 
Principal  Investigator:  Dr.  Norbert  NESSLER 

Thesis  graduates:  Erich  BERLOFFA2^ 

Dieter  FRITZ 
Richard  KLEIN2 > 

Helmut  LEONHARD 
Michael  OBERLADSTXTTER1 } 
Manfred  WIESER 


'  Started  in  1972 
2^  Finished  his  studies  in  1972. 

Scientific  consultant:  Dr.  E.  SALPETER 
Technical  consultant:  Walter  GRADL 

Office:  Brigitte  NESSLER 


in-  i.xmini  iiirtiinimmM 


IV 


ABSTRACT 


The  results  of  the  research  on  the  production, 
transmission  and  reception  of  VLF  and  LF  waves  is 
represented  in  9  compact  chapters. 

A  new  improved  equivalent  circuit  diagram  for  the 
helix  is  calculated  giving  a  good  explanation  of 
the  measured  effects  over  a  large  frequency  range. 

The  receiving  characteristic  of  a  ferrite  rod  an¬ 
tenna  in  the  elliptically  polarized  field  is  in¬ 
vestigated  specifying  a  method  that  allows  all 
allipse  parameters  to  be  measured  simultaneously. 

The  phase  measurement  at  VLF-wave  propagation  and 
the  device  built  for  this  purpose  are  described  and 
an  improved  universally  usable  digital  phase  measuring 
equipment  for  a  large  frequency  range  is  particularized. 

The  various  phenomena,  e.g.  reflection,  refraction, 
and  wave  guiding,  that  are  occurring  at  VLF  and  LF 
wave  propagation  are  described  and  explained. 

A  detailed  description  of  a  long-time  registration 
of  commercial  VLF  transmitters  is  given. 
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CHAPTER  1 

HELICAL  ANTENNAS  IN  THE  VLF- RANGE 


The  Impedances  of  the  antennas  H  1136 ,  H  328,  and 
H  329  Cull  In  air  showed  the  same  characteristic  shape 
of  real  and  imaginary  part  in  the  graphical  representa¬ 
tion  as  has  been  stated  earlier  with  the  antenna 
H  772  /l/.  The  shape  of  the  curves  computed  at  Innsbruck 
University  reminds  of  the  impedance  characteristics  of 
a  quartz  crystal.  It  was  now  near  at  hand  to  examine  the 
equivalent  circuit  diagram  of  a  quartz  crystal  more 
thoroughly  with  particular  consideration  of  its  appli¬ 
cability  for  the  special  case  and  possibly  to  modify  it. 

The  equivalent  circuit  diagram 


Fig.  1.2 

all  CSi  were  united  to  C;  R,  L,  C,  and  Cp  remain  con 
stant.  Fig.  1.3  shows  the  modified  diagram. 


Fig.  1.3 


This  diagram  is  identical  with  the  equivalent  circuit 
diagram  of  a  quartz  crystal. 

Calculation  of  the  impedance  of  the  circuit 
The  complex  conductance  amounts  to 

1 

Y  =  jwCp  + 


R  +  j  (»L  -  -i;) 

wC 


with  co  =  2tt f 


IdtuMUH  IJlIlMiMMMMlIMliaLMaiM.  -■ - — —MMjMMMttMftMl 


3 


1 


or 


R2  +  («L  -  -i-) 


+  j 


tdC 


u)L  — 


zc 


l  R  +  (“L  -  Sc> 


WCP 


(1.1) 


with  real  and  imaginary  part  separated. 


Therefrom  real  and  imaginary  part  of  the  complex  total 
resistance  Z  is  calculated  to  be 


ReZ 


ImZ 


1+  u2Cp2  (R2  +  H2 )  -  2uCp  *  H 

H  -  uCp (R2  +  H2) 

1+  Cp  (R  +  H  )  -  2  Cp  .  H 


(1.2) 


(1.3) 


with  the  abbreviation 


H  *  (uL  -  ) 

«C 


Resonance  frequencies 

The  resonance  frequencies  are  calculated  from  ImZ  *=  0, 
resulting 

H  -  «Cp (R2  +  H2)  *  0  (1.4) 

and  therefrom  we  obtain  a  biquadratic  equation  for  u: 

1 


4  2 


LCp  L 


r2  „  n 

-T  +  2  .  — 


LC 


1  1  1 
+  C  LzCp  +  L2C2 


(1.5) 


with  the  solution 


u 


reS  2LCr 


p2  .  ,  R4Cd2  4R2  Cp2  2R2Cp 

£—  +  _L  +  — i —  /i  +  — -IL  - E-  - - 

T  2 


2L 


LC  2LCr 


LC 


L 

(1.6) 


With  the  antennas  in  question,  R  and  C  are  small,  L  is 
large;  so  the  following  approximations  can  be  introduced: 

2R2Cd  4R2Cp2 

- 2  <  i .  - <  i . 

L  LC 

R4Cn 

The  root  is  expanded  into  a  series,  the  term  — can 
be  neglected. 


resulting 


2RzCp2 

LC 


(1.7) 


(1.8) 

(1.9) 


For  the  special  case  R  =  0  all  neglected  terms  disappear 
so  that  the  resonance  frequency  reads  as  follows: 


2  C+Cp  2  1  z  C 

“i  <R=o)  =  -r-f-  =  “io  *  7=  /  1  +  tt 

LCt-p  /LC  (-p 


(1.10) 


(l.U) 


Thus,  two  neighboring  resonance  frequencies  appear: 

The  first  (lowest)  resonance  originates  in  the  series 
circuit  L-R-C  (Fig.  1.3)  with  Cp  as  capacitive  load 
parallel  to  the  generator. 

The  second  (slightly  higher)  resonance  forms  in  the 
parallel  circuit  L-R  on  the  one  hand  and  C  in  series 
with  Cp  on  the  other  hand.  The  tap  occurs  at  the  capa¬ 
citive  divider  C,  Cp. 


From  (1.8)  and  (i.9) 


(1.12) 


(1.13) 


is  calculated. 

R^Cn 

Generally,  —  <<  1,  i.e.  the  shift  of  the  resonance 
2L 

frequency  on  account  of  the  loss  resistance  is  negligible 
with  the  antennas  used  here  (smaller  than  the  inaccuracy 
of  measurement) . 


Loss  factors  are: 
for  parallel  resonance 


tan  6n  = 


IR 


fR2CE 

2L 

R 


<<  1 


r 


Iojp.L 


wPo,L 


r/ 


L(i  +  f» 
Cp 


RZCp 

for  series  resonance  -it-  «  1 


2L 


tan  5 «  =  T 
s  “S‘L 


wSq  * 


fi 


(1.14) 


(1.15) 


With  the  assumption  of  small  loss  factors  the  resonant 
impedances  are  calculated  as 


R'“PocP>2 


Cp'L 


RCp(c+cp) 


=  Pp 


(1.16) 


for  parallel  resonance,  and 


-  6 


Y  *  5  11  *  j"‘oVR>  *  5 

for  series  resonance. 


i*r 


=  BjS  S  (1.17) 


From  the  loss  factors  both  quality  and  band  width  can  be 
calculated : 


Qs  *  tan  «s  =  i  , 


T  R 


b  c  if  =  tarn  6  •  f 


res 


b.  =  b_  =  -  (1.18) 

S  P  2ir  •  L 


The  calculation  of  the  constants  L»  R,  C,  and  Cp  from 
the  geometric  and  electric  data  of  antenna  and  feed  line 

Inductance  L;  calculated  in  approximation  as  inductance 
of  a  long  air-core  coil  /2/.  An  accurate  calculation  is 
given  in  /3/: 

Inductance  of  a  ningle-layer  cylindrical  coil: 

_  21.8*  n2  *r  ,  -7  ...  .  ... 

L0  =  . .  , _  •  10  |H,  mj  (1.19) 

1  +  2.2  •  — 

r 

21  ....  total  length 


Q  P  0  <j>  p 

-*  Q  r  i  s  *■  i  2?i  *  - 


ooooo 


Fig.  1.4 


If  the  clear  distance  a  between  two  neighboring  windings 
is  not  equal  to  the  wire  thickness  2pj  =  d,  the  correc¬ 
tion  term  AL  must  be  added  to  Eq.  (1.19). 

AL  *  2n  •  Ke  •  r  •  10“7  |H,  m|  (1.20) 

The  correction  factor  K8  can  be  taken  from  Fig.  1.5. 

Series  loss  resistance  Rtotal: 

The  series  loss  resistance  R  of  the  antenna  consists  of 
the  following  partial  resistances: 

Rtotal  “  ^ire  +  ^eddy  curt.  +  Rdiel.+  ®rad. 

1*  R^j-e  is  the  Ohmic  resistance  of  the  wire  where  a 
current  I  with  the  frequency  f  flows.  Current  displace¬ 
ment  by  skin  and  proximity  effect  must  be  considered  at 
the  calculation  of 

The  Ohmic  resistance  Rq  of  the  wire  material  can  be  taken 
from  tables  /4/.  The  Ohmic  resistances  for  wire  gauges 
relevant  for  helical  antennas  are  listed  in  Table  1.1. 

Table  1.1:  Resistance  of  a  round  Cu  wire  in  dependence 
of  diameter. 


0  |mm| 

q  |mm2 | 

R(01  ,/m) 

1  per  R 

0.9 

0.6362 

0.0281 

35.6 

1.0 

0.7854 

0.0227 

43.9 

1.1 

0.9503 

0.0188 

53.2 

1.2 

1.131 

0.0158 

63.2 

1.5 

1.767 

0.0101 

99 

The  increase  of  resistance  caused  by  skin  effect  is 
thoroughly  dealt  with  in  /3/ ,  approximation  formulas 
are  given: 


For 


vs 

—  < 
b 

2 

R  * 

*o 

(1.21) 

2  < 

§<< 

R  = 

*oH^n 

(1.22) 

d 

b  > 

4 

R  * 

r°  a + °-25] 

(1.23) 

d 

b  > 

10 

R  * 

(1.24) 

d  ....  wire  diameter 

b  ....  equivalent  conducting  layer  thickness 

(cf.  Eq.  1.25) 


b.4  Ki  _2 

b  *  - -  •  10  2 

/  f 

|m,  Hz| 

(1.25) 

The  correction  factor  Kj 
material  and  can  be  seen 

is  dependent  upon 

from  Table  1.2 

the  wire 

Table  1.2 

Material 

Ki 

silver 

1.00 

copper 

1.03 

aluminum 

1.37 

zinc 

1.97 

brass 

2.2 

platinum 

2.62 

or  can  be  taken  directly  from  Fig.  1.6  in  dependence  on 
the  frequency. 
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The  resistance  of  a  1  in  silver  wire  in  dependence  on  the 
frequency  and  the  radius  of  the  wire  can  be  deduced  from 
Fig.  1.7.  If  we  use  another  wire  the  effective  radius 
must  be  divided  by  the  factor  K1 . 

Proximity  effect:  If  currents  are  flowing  in  neighboring 
wires  in  same  or  opposite  directions*  current  displace¬ 
ment  occurs  in  addition  to  tne  skin  effect  on  account 
of  the  change  of  magnetic  fields  according  to  Fig.  1.8. 

oo  o  o 

jame  current  opposite  current 

direction  direction 

Fig.  1.8:  Proximity  effect 

The  hf-resistance  of  the  respective  wire  is  thus  increased. 
By  multiplying  the  hf-resistance  by  the  correction  factor 
K2  this  increase  is  taken  into  account.  For  the  case  of 
two  parallel  wires  with  circular  cross  section,  K2  can 
be  taken  from  Fig.  1.9. 

Reddy  current  is  caused  by  the  surrounding  conductive 
rock  and  by  the  influence  of  the  measuring  devices. 

If  the  antenna  is  surrounded  by  a  second  shorted  wire 
with  the  resistance  Rz *  Reddy  current  reads 

M2 

Reddy  current  “  “  ^7  (1.26) 

M  ....  mutual  induction 
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for  the  frequency  range  up  to  about  0.5  MHz,  i.e., 
^ddy  current  ls  increasing  proportionally  to  f2. 

For  higher  frequencies,  the  approximation 


^ddy  current  =  R2  *  (1.27) 

k  ....  coupling  coefficient 


is  valid,  i.e.,  the  increase  of  eddy-current  resistance 
As  only  proportional  to  the  increase  in  R2  by  the  skin 
effect. 


^iel*  c°nsists  of  losses  caused  by  the  capacitive 
currents  in  the  wire  isolation  and  in  the  carrier  ma¬ 
terial.  For  constant  current  R^^  increases  with  f 3 . 

Rdiel.  *  0.25  •  f3  •  CL  •  L2  •  dc  .  lo3  (1.2 

|n»  HZ,  F,  Hj 

CL  ••••  self-capacitance  of  the  antenna 

L  .  inductivity  of  the  antenna 

dc  ....  loss  factor  of  the  dielectric 

The  value  of  CL  can  be  seen  from  Fig.  1.10  (with  the 
same  designation  as  in  Fig.  1 . 4)  ,  dc  froir  tables, 
e.g.  in  /4/. 


^radiation  must  be  considered  only  with  very  high  fre¬ 
quencies  or  in  the  case  of  antenna  dimensions  larger 
than  A/lo  (A  =  working  wavelength),  in  /5/  radiation  re 
sistances  for  smaller  antennas  (capacitive  or  inductive 
standard  antenna)  are  quoted: 


E 

o  tc  E 

*-  XL 

O 


o 


I 

o 


N 

t 


L'l  Old 
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capacitive 

antenna 


Rg  *  80  •  u2 


inductive 

antenna 


R*  -  80  •  »2 


(1.29) 

(1.30) 


X  .  wavelength 

1  .  length  of  the  electric  dipole 


n«A  ...  turn  area  of  the  magnetic  dipole 
«e'*m  .  radiation  resistances 


In  every  case,  however,  R^d  is  essentially  (several 
decimal  powers)  smaller  than  Ryire  and  Rdiel'  respectively. 

With  sufficiently  large  gallery  cavity,  Reddy  current  is 
as  well  much  smaller  than  and  Rdiel  on  account  of 

the  low  conductivity  of  the  surrounding  rock  and  the 
small  dimensions  of  the  used  measuring  devices. 


Capacity  C:  For  the  calculation  of  C,  the  approximation 
formula  for  the  static  capacity  of  a  symmetric  cylindrical 
antenna  has  proved  useful.  According  to  /6/ 


3.54  •  1 


In 

i ' 

-  -  0.55 

10~13  | F,mj 
(1.31) 


21  ....  total  length  of  the  antenna 

2 1  >  r 

r  .  radius  of  the  cylinder 

In  this  simple  form  the  equation  is  only  valid  with 
comparatively  narrow  windings  (as  e.g.  at  H  772,  H  1362, 
H  1136  Cull,  H  328  etc.) . 


Cp  is  calculated  from  the  parallel  connection  of  the 
capacity  of  the  feed  line  and  the  capacity  of  the  two 
inner  windings: 


Cp  =  CK  +  Ci 
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According  to  /3/  the  feed-line  capacity  appears  as  the 
capacity  of  the  symmetric  two-wire  line  and  reads 


0.121  .  er  .  1 


.10  ,  , 

10  |F,m| 

(1.32) 


D  ....  distance  between  wires 
d  ....  diameter  of  the  wires 
1  ....  length  of  the  wires 


Generally,  the  cable  capacity  is  given  by  the  manufacturers 
or  can  easily  be  determined  by  experiment. 

The  capacity  of  the  two  inner  windings  of  the  antenna 
can  be  calculated  according  to  /3/  as  the  capacity  of 
two  parallel  wires  of  the  length  1,  of  the  antenna 
circumference  2r* ,  the  distance  D  and  the  diameter  d: 


0.12l*c.»l 


0.754*er»r 


lo"10  |F,m| 


(1.33) 


with  er  to  be  determined  experimentally  on  account  of  the 
influence  of  t* ±  wire  isolation  U  not  known  with  suffi¬ 
cient  accuracy) . 
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quantity 

unit 

H  1136  Cull 

H  328  I 

H  328  II 

Remark 

n 

1136 

328 

328 

r 

m 

1.43*10_1 

5.5*10"2 

5.5*10"2 

8 

m 

1.72* 10”3 

2.5«10"3 

2 . 5* 10”3 

p 

m 

0.75* 10”3 

5*10"4 

•4 

5*10 

* 

1 

m 

0.977 

0.41 

0.41 

loss  factor 

*c 

4 • 10” 3 

4* 10” 3 

n 

1 

O 

H 

• 

of  the  dielectric 

carrier  material 


Table  1.4s  Results 


quantity 

unit 

equation 

H  1136  Cull 

H  328  I 

H  328  II 

L 

H 

19,  20 

2 .7  *10”2 

8  .04  •  ICf4 

8. 04  •lo”'* 

Kwire 

n 

24,  25 

57 

16 

16 

+  *diel. 

a 

28 

332 

24 

24 

*  Rtotal 

n 

389 

40 

40 

C 

F 

31 

20.1*10"12 

7.8-10"12 

7.8-10"12 

F 

32,  33 

90*  10"12 

98- 10"12 

1.97*10"9 

H  328  II;  Same  data  as  for  H  328  I,  cp  artificially  in¬ 
creased  by  parallelling  of  a  low-induction 
ceramic  drop  capacitor. 
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The  measurements  were  executed  in  a  larger  cavity  in 
the  SO-gallery  of  the  mine  St.  Gertraudi/Tyrol .  The 
following  measuring  devices  were  used: 

for  the  frequency  range  up  to  600  kHz: 

selective  levelmeter  SPM-2,  Wandel  &  Goltermann 
decade  R-C -generator  TG  66  A,  Levell 
three-stage  preamplifier,  gain  80  dB,  50-1000  kHz, 
transistorized,  constructed  in  the  VLF-lab 

for  the  frequency  range  over  600  kHz: 

L-C-generator ,  transistorized,  constructed  in  the 
VLF-lab 

selective  heterodyne  voltmeter  2006,  Brtiel  &  Kjaer 

for  both  frequency  ranges: 

R-C-L-bridge  LE  300/A  1,  Hatfield 
frequency  counter,  type  100A,  Monsanto 


The  values  of  RE  Z  (exp)  and  IM  Z  (exp)  were  computed 
according  to  the  formulas  given  in  /l/  and  listed  in 
tables  1.5  to  1.7. 

The  columns  mean  from  the  left  to  the  right: 


f 

kHz 

working  frequency 

RE  Z  (exp) 

ka 

real  part  of  impedance  calculated 
from  experiment 

RE  Z  (th) 

k« 

real  part  of  impedance  calculated 
from  the  equivalent  circuit  diagram 

DRE  Z 

ratio  RE  Z  (exp)  /  RE  Z  (th) 

IM  Z  (exp) 

k  a 

imaginary  part  of  impedance  calculated 
from  experiment 

IM  Z  (th) 

kn 

imaginary  part  of  impedance  calculated 
from  the  equivalent  circuit  diagram 

DIM  Z 

ratio  IM  Z  (exp)  /  IM  Z  (th) 

A  graphical  representation  of  the  measuring  results  is 
shown  In  Figs.  1.11  to  1.16. 


r*  ®  IT  IT,  Xt  ®  O'  (VI  « 

O'  C5-  .:  J  _  o-  Lt.’  ir  01  ®  ®  r- 

X  »•••«••••*  ~ 

M  COOOOOCOOC;'  >000000 

o 


X) 

<0 

I- 


(\i  (Vi  (Vj  o  o  o  r\j  m 


X 

O 

0  0 

O  O 

O 

O 

0 

O 

000 

O 

00000 

♦ 

♦  4 

♦  ♦ 

♦ 

♦ 

♦ 

♦ 

♦  ♦  ♦ 

♦ 

♦  ♦  ♦  ♦  ♦ 

■w 

U- 

U.  UI 

UiUJ 

U J 

OJ 

UJ 

UJ 

UJ  UJ  UI 

UI 

Lil  Ui  111  Ui  UI 

4  ® 

«in 

p-H 

0 

(Vi 

-h  in  in 

0 

HHOrtrt 

® 

(V  r~ 

®  ■£) 

r- 

0 

in 

<£i 

r^MO 

® 

IM 

i*i  0 

O  O' 

a 

•0 

IV 

<-H  O  ® 

(V 

4  <0  O'  4  K 

«r 

<v< 

««•  f— « 

« in 

r> 

IV 

H 

® 

(V  ®  — ' 

(Vi 

♦  ♦nnm 

M 

0 

O  O 

OO 

O 

0 

O 

c 

OOO 

0 

00000 

1 

l  l 

t  • 

• 

« 

t 

i 

( 

1 

1  1  »  1  t 

(VI 

(V  (Vi 

r-*  «mi 

p*4 

p-4 

r-4 

O 

f"4 

fv» 

<v 

(H  ®-*  ^ 

CL 

O 

O  O 

O  O 

O 

O 

O 

O 

O 

O 

O 

O 

OOOOO 

X 

♦ 

♦  4 

♦  4 

♦ 

♦ 

♦ 

♦ 

♦ 

1 

♦ 

♦ 

♦  ♦  ♦  ♦  ♦ 

UI 

UJ 

Ui  UI 

UI  UI 

UJ 

UJ 

UJ 

Uj 

UJ 

Ui 

UJ 

Ui 

UJ  Ui  Ui  Ui  IaJ 

m 

tt  IT 

in  ® 

0 

0 

® 

0 

0 

c 

<0 

4 

amioif.  w 

o>  n 

0  ® 

o> 

CO 

#0 

fO 

® 

ro 

•»4 

nnroHfi 

IM 

® 

m 

4  0 

iH 

® 

a 

(VI 

O' 

in 

r> 

(MoaniK 

(VI 

«  ® 

* 

(VI 

r-4 

r- 

CO 

(VI 

n 

r-  in  4  4  i*» 

X 

• 

•  • 

•  • 

• 

• 

• 

• 

n 

• 

• 

0 

O  O 

OO 

O 

0 

O 

O 

O 

CJ 

O  O 

0 

0 

0 

0 

0 

1 

l  l 

»  1 

( 

1 

1 

l 

1 

1 

1 

»  1  »  1  1 

M 

0 

cn  ® 

(VI  (V 

O 

sr 

4  n  (V  ®  0  Ch 

K  D  O  4) 

iH 

n  a* 

r*-  *-1 

<0 

00® 

*-HI 

(vi-io^irfi^n 

UI 

• 

•  • 

•  • 

• 

• 

• 

• 

• 

• 

• 

• 

U 

4 

•-4  •— « 

0  ♦ 

f—4 

.-■4 

O 

O 

«~4 

^•4 

om  (vi  ®  in 

c 

(V 

(*)  (Vi 

r-4 

x  ♦  id  x 

p-4 

—• 

— • 

O  O 

0  0 

0 

0 

0 

0 

0 

0 

r-4 

*•4 

^  O 

O 

O 

O 

X 

O 

O  O 

O  O 

O 

0 

O 

0 

0 

O 

O 

O 

O  O 

O 

O 

0 

»- 

♦ 

♦  ♦ 

♦  4 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

1  * 

♦ 

♦ 

♦ 

*- 

Uj 

Ui  u 

Ui  UI 

UJ 

Ui 

UI 

UI 

Ui 

Ui 

UI 

uJ 

UI  Ui 

lu 

Ui 

Ui 

® 

4  cm 

®  (VI 

in 

f- 

4 

r> 

0 

m 

CM 

*  0 

p- 

♦ 

O' 

4 

N  « 

KI9 

(*> 

(VI 

4- 

(VI 

0 

O 

in  (*- 

tn 

<0 

« 

IM 

0 

in  ® 

n  m 

® 

(VI 

® 

O 

(VI 

O'  c- 

CO 

♦ 

m 

m  it 

rxvi 

1—4 

H 

(VI 

(VI 

(*> 

m 

0 

0 

r-4 

r-4 

*■4 

a 

O 

0  0 

O  O 

0 

O 

0 

0 

0 

0 

O 

O 

O  O 

O 

0 

0 

*X« 

^4 

^4 

0 

O 

O 

0 

0 

0 

^4 

—4 

(Vi 

(Vi 

(VI 

a 

O 

0 

0 

O 

O 

0 

0 

0 

0 

0 

0 

O 

0 

0 

0 

X 

> 

» 

1 

♦ 

♦ 

4 

♦ 

♦ 

♦ 

♦ 

♦ 

» 

1 

» 

1 

t 

UI 

UI 

U; 

•j 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UI 

UJ 

UI 

UI 

Ui 

fw 

•* 

4 

4) 

*•4 

in 

0 

0 

r- 

in 

a> 

0 

0 

4 

(VI 

in 

V 

<x> 

* 

n 

m 

0 

0 

* 

4- 

O 

IM 

m 

r-4 

r*-% 

n 

X 

n 

<* 

r- 

r-4 

CO 

<-* 

(Vi 

H 

(VJ 

n 

u 

*■« 

(V 

r> 

0 

flO 

in 

■+ 

(VI 

H 

UI 

• 

• 

• 

■ 

• 

> 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

a 

O 

O 

0 

0 

0 

0 

') 

0 

0 

0 

0 

0 

O 

0 

0 

O 

0 

0 

0 

0  0 

O 

O 

0 

O 

0  0 

,0 

,0 

.0 

.0 

O 

0 

0 

0 

<r  0 

O 

O 

If,  or. 

V>  (V 

0  0 

0000 

O 

•S' 

0 

(VI 

m  ® 

O 

p"4 

p"4 

(VI 

^  sr 

0  m  0  »n 

0 

u 

pH 

r-H  p-H 

(VI 

C\J 

(VI  .1 

•\i  i  n 

n  n  4  4 

in 

>1 


a 


[Mb 


j 

j 

i 


i 


/</-  b 


N  MC(V.Cj|04MT(UI>irrw^Al^4(0 

hcc-io(v.r-ar>o^nO'inuiM 

2  . . . 

IV  IV  !\J  r*  ^  *H  — »  ri  C  O  f**“  (V  IV  IV 


I  O  O  O  «-!  •-• 

>000000 
♦  ♦  4  ♦  I  I 
I  lu  U  ll<  UJ.Uj  UJ 

1  tv  tv  iv  -#  r--  in 

1  (D  >4  V  4  ID  1*1 
inniHtfi  4  4 
1  ®  ♦  m  s  <*) 

>•••••• 

>000000 

itiii 


ooohuih, 
o  o  o  o  o  o  < 
****** 
uj  uj  ul  ui  ui  taJ  1 
0404 
r*  IV  *“•  8)  8)  I 

®  IV  O'  ®  »  ■ 

1*1  ®  IV  H  1 
•  ••••• 
OOOOOOI 
I  I 


->ocooo<-><^co^-»4>hooo 
eeoooooooooooooo 
******  1  !♦♦♦♦♦♦♦♦ 
UJUJuiUJUIIUUJUIUJUJUJUJkJUIUJUJ 
D«M*l^4aa)4SCM)*ON 
«C(\IOO«ff'OMJ404NON 

4HITH00404  o  tn  <C  CO 

OC  OOOOOO  OOOOOOOO 

I  I  I  I  I  1  I  I  I  I  I 


•*>  ~*  «C  »  ®  O'  OIM  OlP  «  «M\|N 
(fNOO»0'»OB4>rtO(M\;Nr> 
•  ••••••••••••••a 

m  oooooo^hoooo-<o^p~ 

O'  r- 

<v 


IVI*)I*)4  H>a*<»<i»l>aH>ar<(,)r«l|| 
OOOOOOOOOOOOOOOO 
I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 

UiUiUJUJUJUJUJUJUJUJIiJUJUJUJUJUi 

nu(j'4M?njDn(vni:in»->4 

tn  04  £  r-.  ir.  -»r^®m^inh-0* 
inoo>NH4  £  £  oon<x>®o~*® 
r*  h  o  »h  (v  n  ^inininivinmiva 

OOOOOOOOOOOOOOOO 


oi,tMyj'*»<«"««o«*<ttii,i'4 
COCOCOOOOOOOO  00  o 
1  I  I  I  I  1  I  I  I  I  ♦  ♦  I  I  I  I 
UJUjUJUjUJUJUjUJUjUJUJUJUJUJUJUJ 
IM  (\m  l\l  4  «  CD  4  8)  O  81  4  (VIN 

a.  '.^c\jr»£®r^n(v^®«ino>® 

o  C  5  &  *  «0fl^4flB»4«nn 

&  t  4.'  >■+  fv  n  ■*  •c  ®  1*1  i*j  iv  iv  f» 


0000 


COOOOOOOOOO 


OOOOOO  OOOOOOOO 

•  •••••«  •••••••• 

OOOOOOOC-  OOOOOOOO 

occomoMn^idiiocooo 

Clf.rtJffOOCOOOHMnoO 


15 


Discussion  of  results 

The  deviations  still  appearing  between  experiment  and 
theory  can  be  explained  on  account  of  the  following  con¬ 
siderations: 

1.  The  equivalent  circuit  diagram  is  calculated  with  a 
constant  resistance  Rtotal  which,  however,  changes  with 
frequency  (skin-effect,  dielectric  losses) .  The  influence 
of  the  surrounding  medium  on  the  eddy-current  resistance 
is  virtually  incalculable. 

2.  Due  to  the  nigh  humidity  in  the  mine  the  parameters 
Rdiel  and  C  change  in  a  way  not  yet  investigated  (measure¬ 
ments  in  the  mine  have,  however,  among  others  the  advantage 
of  an  essentially  lower  noise  level! 

3.  The  parameters  L,  R,  C,  Cp  cannot  be  calculated  with 
the  required  accuracy  (approximation  formulas)  because 
there  are  quantities  (e.g.  t)  that  are  not  known  exactly 
enough. 

4.  The  calculation  of  RE  Z  (exp),  IN  Z  (exp)  is  subject 
to  an  error  of  measurement,  yet  not  exceeding  10  %. 

The  equivalent  circuit  diagram  gives,  however,  a  good 
interpret.-!  Ion  of  some  effects,  that  have  not  been  able 
to  be  clearer  tisfactorily  up  to  now: 

1.  In  the  equations  of  /l/  Cj.  is  contained  in  the  cal¬ 
culation  of  the  series  ..  ^sonance  frequency.  This  is,  how¬ 
ever,  in  variance  with  present  experiments,  where  the 
decrease  of  Cj_  by  raising  th.  c.  stance  between  the  two 
inner  windings  has  no  influenc-..  i  the  series  resonance. 

The  equivalent  circuit  diagram  shows  definitely  that 
is  only  affecting  the  parallel  resonance.  This  is  proved 
by  the  experiment  with  the  antenna  H  328  II,  with 


16 


being  artificially  increased  by  parallelling  a  capacitor. 
The  series  resonance  frequency  remained  constant  whereas 
the  parallel  resonance  approached  the  series  resonance. 

2.  Body  effect:  The  small  static  capacity  C  is  essan- 
tially  changed  (increase  of  efl)  by  approaching  a  hand 
especially  at  the  two  ends  of  the  antenna,  whereas  it  is 
relatively  insensitive  in  the  middle. 

3.  Influence  of  humidity:  Water  vapour  deposits  on  the 
antenna  on  account  of  the  high  humidity  thus  influencing 
both  loss  resistance  and  capacity  C.  This  gives  a  good 
explanation  for  the  difficulties  occurring  with  earlier 
experiments  using  wood  as  carrier  material.  Water-re¬ 
pellent  materials  such  as  epoxy,  PVC  etc.  proved  more 
favorable. 

Finally  the  qualitative  agreement  of  the  impedance  curve 
in  experiment  and  theory  can  be  regarded  as  best  proof 
for  the  validity  of  the  equivalent  circuit  diagram  which, 
however,  holds  only  in  the  range  of  the  first  (deepest) 
series  and  parallel  resonances.  At  higher  frequencies 
additional  effects  of  the  individual  windings  occur  and 
the  simple  model  described  here  fails. 

Considerations  on  the  efficiency  of  the  antennas  in  air 

According  to  /!/ ,  the  distant  field  of  a  small  self¬ 
resonant  helix  with  many  windings  may  be  regarded  as  a 
combination  of  an  electric  with  a  magnetic  dipole.  The 
length  of  the  helix  corresponds  to  the  length  of  the 
electric  dipole,  the  turn  area  corresponds  to  the 
magnetic  dipole  n»A  =  n*r2,n. 

For  these  dipoles  the  efficiency  can  be  calculated  with 
given  radiation  frequency  (here  the  first  resonance  fre¬ 
quency)  according  to  /5/: 


radiation 


radiation 


Prad  +  ploss  l2 ^Rrad  +  Rioss^ 


(1.34) 


As  mentioned  already  the  radiation  resistance  is  smaller 
than  the  loss  resistance  by  several  decimal  periods,  so 
that  Rrad  can  be  neglected  in  the  denominator: 

l2prad  Rrad  prad  ,,  ,ri 

n  =  -  = - —  =  --—■■■■  (1.35) 

I2r1oss  Rloss  r1oss 


Rrad  can  be  calculated  according  to  Eqs.  (1.29)  and  (1.30) 
Eqs.  (1.34)  and  (1.35)  are  only  valid  with  antenna  dimen¬ 
sions  «  X  (which  is  certainly  the  case  for  the  antennas 
described).  Table  1.7  gives  a  survey  of  these  calcula¬ 
tions  for  the  antennas  H  1136  Cull,  H  328,  and  H  1136 
in  water. 

Table  1.7 


Rloss 

n 

Re 

n 

dB 

n 

dB 

fres 

kHz 

H  1136  Cull 

400 

1 . 6  *  10~3 

-74 

40-10~6 

-70 

220 

H  328 

40 

0.24 

-22 

5.6*10~4 

-48 

2020 

H  1136  Cull 
in  water 

35 

00 

1 

o 

r-4 

• 

• 

-71 

2.8*10'6 

-89 

lcl 

Discussion:  Antennas  with  a  relatively  low  resonance  fre¬ 
quency  in  air  (e.g.  h  1136,  H  772)  show  a  decrease  of 
resonance  frequency  caused  by  a  narrower  winding  and 
thus  increased  inductance.  As  regards  efficiency  these 
"inductive"  antennas  are  by  far  less  favorable  than 
"capacitive"  antennas  as  is  shown  in  the  table.  Thus, 
for  any  application  we  have  to  make  a  compromise  between 
antenna  efficiency  and  favorable  frequency  for  the  trans¬ 
mission  in  the  dissipative  medium,  which  is  thoroughly 
dealt  with  in  /8 /. 
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Diagram  1.17  shows  the  characteristic  course  of  a  propaga¬ 
tion  measurement.  In  this  case,  antenna  H  1136  Cull  was 
operated  in  water  in  the  mine  St.  Gertraudi. 

Curve  a  shows  the  curvature  of  field  strength  in  dependence 
on  the  distance.  The  antenna  was  situated  in  the  water- 
filled  "Mabuse-gallery" ,  measurements  were  executed  along 
the  west-gallery. 


Measuring  data  of  this  experiment: 
frequency  f  -  1.1  kHz 

antenna  current  I  -  1  A 

applied  power  =  100  W 

Transmission  was  detected  unobjectably  in  a  distance  of 
1  km  with  a  signal-to-noise-ratio  of  7  dB.  The  com¬ 
ponents  Br ,  could  be  determined.  BQ  vanished  in  the 
noise  at  greater  distances. 

This  result  indicates  the  predomination  of  the  part  of 
the  electric  dipole. 


Stimulated  by  this  experiment  the  same  antenna  was  once 
more  fed  in  water,  this  time,  however,  in  a  lower  water- 
filled  gallery  near  the  VLF- laboratory.  For  this  experi¬ 
ment,  the  VLF -transmitter  (Savage)  in  the  lab  was  used. 

Measuring  data  of  the  experiment: 
frequency  f  =  1.1  kHz 

antenna  current  I  =  3 , 5  A 

applied  power  Pjl  =  700  W 

The  power  could  not  be  increased  on  account  of  sparkovers 
at  the  ends  of  the  antenna  that  would  surely  have  in¬ 
fluenced  the  measurement. 

Surprisingly,  the  signal  was  to  be  detected  unobjectably 
even  in  a  distance  of  more  than  7  km  on  the  earth  surface. 


CHAPTER  2 


MEASURING  THE  ROTARY  FIELD  ELLIPSE  BY  A 
FERRITE  ROD  ANTENNA  -  "PEDAL  LOCUS  CURVE" 


In  the  magnetic  alternating  field  a  coil  (frame)  receives 
an  e.m.f.  being  proportional  to  the  area  A,  the  number  of 
turns  n  and  the  change  velocity  of  the  magnetic  induc¬ 
tion  B,  the  latter  being  proportional  to  the  projection 
of  the  field  direction  to  the  frame  axis.  The  desired 
high  e.m.f.  is  reached  by  a  high  p  of  the  ferrite  mate¬ 
rial  ,  filled  into  the  coil  and  A.  Technical  problems  re¬ 
strict  the  two  values.  Comprehensive  investigations  /9/ 
showed  the  use  of  a  ferrite  rod  to  be  optimal  with  the 
dimensions  given.  This  ferrite  rod  coaxial  to  the  re¬ 
ceiving  coil  represents  the  coil  axis  and  thereby  the 
direction  of  the  receiving  antenna  for  all  further  con¬ 
siderations. 

An  alternating  field  constant  in  space  induces  an  e.m.f. 
in  the  ferrite  antenna  according  to 

e.m. f . in(jt  (t)  =  K*B  sinwt  cost (2.1) 

K  ...  calibration  factor,  specific  for  electric  and 
geometric  properties  of  the  ferrite  antenna 
4*  ...  angle  between  ferrite  rod  and  field  direction 

If  the  field  consists  not  only  of  one  vector,  but  is  an 
elliptical  field,  then  the  projection  of  the  instantaneous 
field  vector  on  the  antenna  direction  is  a  direct  measure 
for  the  e.m.f. 
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2.1  Two-dimensional  case 

At  first,  the  three-dimensional  problem  is  confined  to 
a  two-dimensional  one  by  moving  the  antenna  only  in  the 
plane  of  the  rotary  field  ellipse.  Let  the  x-y  plane  of 
the  coordinate  system  coincide  with  this  plane. 

It  could  be  shown  in  / 11/  that  always  only  one  plana/ 
ellipse  will  result,  even  if  more  than  two  components 
with  different  phase  and  direction  in  space  are  super¬ 
posed. 

On  account  of  the  continuous  change  of  direction  and 
value  of  the  field  vector  the  instantaneous  vector  in  the 
antenna  direction  generally  does  not  yield  the  maximum 
amount  of  projection  and  thus  maximum  e.m.f.  In  Pig.  2.1a 
the  maximum  e.m.f.  is  plotted  as  radius  for  each  position 
of  angle  of  the  receiving  antenna  (dashed  curve  "pedal 
locus  curve")  and  reduced  to  the  same  scale  the  rotary 
field  ellipse  is  plotted.  It  can  be  realized  that  always 
that  point  T  of  the  ellipse  produces  one  point  of  the 
pedal  locus  curve  F  in  which  a  normal  to  the  antenna  di¬ 
rection  contacts  the  ellipse. 

In  vector  form  the  analytic  description  o 2  this  pedal 
locus  curve  measured  in  earlier  studies  /10,11/  reads  as 
follows : 

+  -  COSl|< 

Xp  =  va,cos2i|>+  b2sin2ij/  (2.2) 

simp 

being  the  basis  for  further  calculations. 


Special  cases: 


For  a  =  b  (circular  polarization) 


a 


COSiJj 

sin<J< 
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or  without  parameter 

x2  +  y2  ~  a2  (circle  is  counted  twice) 

The  receiving  e.m.f.  is  independent  on  the  direction 
(see  Pig.  2.1b). 


For  b  =  0  (linear  polarization) 

Z  j  fcos* 

*F  =  a  cos *  [sin*] 


or  without  parameter 


and 


x2  +  y2  ±  ax  =  0 
a  _2 

(x  i  — ) 2  +  y2  =  —  (2  circles  shifted  by  ±a/2 

2  4 

including  the  origin) 


in  the  standard  form,  resp.,  reflecting  the  figure-eight 
pattern  of  the  pure  cosine-dependency  (see  Fig.  2.1c). 


2.2  Phase  relation 


Assuming  linear  polarization  (only  one  vector  instead  of 
the  rotary  field  ellipse)  it  becomes  evident  that  the  in¬ 
duced  voltage  Uj.nd  in  Phase  with  the  alternating  field, 
independently  of  the  angle  1 1>  between  field  vector  and 
antenna  direction  (cf.  Fig.  2.1c).  If  there  exists  a  rotary 
field  ellipse,  each  rotation  of  the  instantaneous  vector 
at  the  same  time  means  a  continuous  phase  rotation  on  2v . 
With  the  phase  position  of  the  major  axis  of  ellipse  as 
zero  point  we  measure  a  phase  for  each  antenna  direction 
corresponding  precisely  to  that  instantaneous  vector  that 
induces  maximum  amplitude  into  the  antenna,.  With  the  excep¬ 
tion  of  the  case  when  the  antenna  is  oriented  in  the  main 
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axes  of  the  ellipse  the  phase  of  the  e.m.f.  in  the  antenna 
is  not  identical  with  the  phase  of  the  instantaneous  vcc-tr. 
in  antenna  direction. 

Corresponding  to  the  definition  of  the  pedal  locus  curve, 
the  phase-determining  contact  points  T  (cf.  Fig.  2.1a)  of 
the  tangent  lines  normal  to  the  antenna  move  close  toget.hc 
near  the  major  axis  but  move  away  from  the  minor  axis. 

This  means  that  the  phase  of  a  dominates  in  quite  a  large 
angular  sector  around  A.  In  the  neighborhood  of  B,  how¬ 
ever,  the  phase  position  changes  much  more  than  in  A  even 
at  the  same  change  of  space  angle.  This  is  of  great  im¬ 
portance  for  the  phase  measurement  described  later  on:  the 
phase  errors  caused  by  errors  of  angular  adjustment  are 
lowest  at  maximum  amplitude. 


2.3  Technology  of  measurement 

For  measuring  the  rotary  field  ellipse,  a  ferrite  antenna 
with  Cardanic  suspension  xs  used.  Fig.  2.2  shows  the  gonio 
meter  that  had  been  used  earlier  and  is  described  in  /!C/. 
This  special  construction  allows  to  turn  the  top-most 
plane  C  (on  which  the  'errite  rod  is  pivoted)  into  any 
direction.  The  scale  has  a  reading  accuracy  of  1/2  degree. 

If  the  position  of  the  ellipse  plane  is  known,  e.g.  let 
it  be  horizontal,  a  and  b  are  comparatively  easy  to  measur 
The  top-most  plane  of  the  goniometer  is  adjusted  in  the 
direction  of  the  ellipse  plane  and  the  minimum  and  maxima 
voltages  are  measured.  For  determining  the  directions  of 
a  and  b,  "minimum"  adjustment  (=  b)  is  preferred,  yielding 
higher  accuracy  (cf.  Fig.  2. la). 

It  is  much  more  difficult  to  measure  the  rotary  field 
ellipse  if  the  latter  is  oriented  arbitrarily  in  space, 
orientation  being  unknown. 


There  are  five  independent,  unknown  quantities  to  be 
determined: 

Orientation  of  ellipse  plane  in  space  2  angles 

Orientation  of  the  major  (or  minor) 

axis  in  the  ellipse  plane  1  angle 

Size  of  axes  2  voltages 

In  /10/ ,  two  different  methods  were  developed  for 
measuring  these  quantities  by  means  of  the  goniometer. 

Assumption  for  method  1:  When  rotating  the  ferrite  rod 
in  a  plane  arbitrarily  inclined  to  the  ellipse  the  mi¬ 
nimum  reception  voltage  always  occurs  in  the  plane  normal 
to  a,  i.e.#  the  plane  given  by  the  z-axis  and  the  minr>r 
axis  b.  Proceeding  on  this  assumption  the  5  desired 
quantities  can  be  measured  as  follows  in  5  definite 
steps  (cf.  Fig.  2.2). 

1.  Starting  position:  plane  A  horizontal,  plane  B  in 

northern  direction  — >  a  =  0°,  y  =  0°, 
ferrite  rod  in  plane  B 

2.  3  =  0°,  y  =  90° ,  turn  a  as  far  as  U^nd  =  min;  fix  a 

3.  y  =  0°,  turn  &  as  far  as  U^nd  =  min;  fix  3_ 

4.  Turn  y  as  far  as  -  min  — ■>  U3 ,  y 

The  ang1  .  ",  y  resulting  therefrom  are  listed. 

5.  y  ±  90°,  -Jj.nd  mm*>  v; 

6.  y  =  0O,  3  ±  90°,  Uind  — >  Ui_ 

Uj  ...  proportional  to  major  axis  a 
U2  ...  proportional  r.o  minor  axis  b 
U3  ...  should  be  zero  u  oise  level) 

A  correction  of  the  angles  a  oi  b  for  an  accurate  adjust¬ 
ment  of  U3  can  be  done  after  step  4.  The  minimum  voltage 
Uo  must  be  as  small  as  possible  compared  to  Jj  thereby  con¬ 
stituting  a  measure  for  the  accuracy  of  the  angular  measure¬ 
ment. 
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As  is  shown  in  the  following  calculation,  the  basic 
assumption  for  method  1  only  yields  an  approximation 
that  is  fulfilled  the  better  the  narrower  the  allipse, 
i.e.  the  higher  the  rat:o  a/b  and  the  smaller  the  incli¬ 
nation  of  the  ellipse  plane  to  the  horizontal  one  from 
which  the  measurement  starts. 

Method  2  is  based  on  the  fact  that  there  is  always  one 
space  direction  where  an  absolute  minimum  of  reception 
voltage  can  be  Measured:  the  normal  to  the  ellipse 
(z-direction)  which  is  reached  within  several  steps  of 
asymptotic  adjustment.  If  the  normal  to  the  ellipse  and 
thus  the  position  of  the  ellipse  plane  is  determined, 
subsequent  measurements  are  performed  as  described  earlier. 

Method  2  is  always  successive  except  for  b  =  0  and  will 
be  preferred  in  such  cases  where  method  1  delivers  in¬ 
accurate  results. 

With  method  1  being  technologically  much  more  favorable 
(5  definite  steps  of  adjustment  on  the  one  hand,  asymp¬ 
totic  finding  of  the  first  adjustment  on  the  other  hand) , 
an  antenna  pattern  is  first  measured  in  the  horizontal 
plane.  Assuming  random  position  in  space  of  the  rotary 
field  ellipse,  the  horizontal  pattern  means  an  inclined, 
plane,  central  sectic,..  through  the  pedaloid. 

For  the  mathematical  description  of  this  measuring 
technology  a  detailed  calculation  of  receiving  charac¬ 
teristics  is  necessary.  In  the  general  case,  the  basic 
element  of  the  curves  and  surfaces  listed  in  the  fol¬ 
lowing  is  the  ellipse. 
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2.4  Pedal  locus  curve  of  the  ellipse 


General  definition  for  a  pedal  locus  curve  (cf.  /12/) : 

The  pedal  locus  curve  of  a  given  curve  related  to  a 
pole  0  is  the  geometrical  locus  of  the  pedal  loci  F 
of  those  perpendicular  lines  that  can  be  dropped  from 
C  on  to  the  tangents  of  the  curve. 

Examples : 

a)  cardioid  -  pedal  locus  curve  of  the  circle  for  one  of 
its  points  on  the  periphery  as  a  pole  (Fig.  2.3a) 

b)  lemniscate  or  figure-of-eight  curve  -  pedal  locus 
curve  of  an  equilateral  hyperbola  with  the  center 
as  a  pole  (Fig.  2.3b). 

For  the  ellipse  the  pedal  locus  curve  takes  the  form  as 
is  illustrated  in  Fig.  2.3c  with  the  center  as  a  pole. 

Two  essential  features  are 

1.  that  the  pedal  locus  curve  differs  the  less  from  its 
ellipse  the  mere  circular  the  latter  becomes,  and 

2.  that  the  pedal  locus  curve  contacts  the  ellipse  in 
the  end  points  of  the  major  and  minor  axes. 

From  the  parametric  representation  in  Eq.  (2.2)  the 
arbitrary  radius  ..  of  the  pedal  locus  curve  can  be 

r 

derived.  This  calculus  was  done  by  /13 /  and  /14/. 

rn  =  /  a2 cos2  a  +  b2sin2  a  =  /  — — cos  2  a 
P  2  2 

(2.3) 

Therefrom,  the  parameter-free  representation  of  the 
pedal  locus  curve  (Fig.  2.3c)  in  Cartesian  coordinates 
can  be  derived: 

a2x2  +  b2y2  =  (x2  +  y2)2  (2.4) 
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2,5  The  trlaxlal  pedaloid 

constitutes  the  three-dimensional  counterpart  to  the 
pedal  locus  curve  of  the  ellipse  and  can  be  derived  from 
the  ellipsoid  on  account  of  the  definition  of  the  pedal 
locus  curve. 

The  derivation  of  the  parameter-free  representation  of 
the  pedaloid  is  done  in  the  same  way  as  with  the 
ellipsoid,  i.e.,  in  the  arbitrary  radius  of  the  pedaloid 

rped  *  ^rp  cos2 3  +  c2  sir.2 8  (2.5) 

rp  is  substituted  by  Eq.  (2.3)  and  the  angular  functions 
by  the  relations  mentioned  above.  The  following  equation 
results: 

a2x2  +  b2y2  +  c2z2  *  (xI:  +  y2  +  z2)2  (2.6) 


Special  case  1: 

revolution  pedaloid  -  it  develops  by  rotating  a  pedal 
locus  curve  of  an  allipse  e.g.  on  its  major  axis.  Since 
b  -  c  in  this  case, 

a2xz  +  b2(y2  +  z2)  =  (x2+y2  +  z2)2  (2.7) 

is  valid. 

Note:  Pedaloid  and  revolution  pedaloid  occur  in  optics 
at  the  description  of  the  double  refraction  of  the  light. 
The  ray  surface  or  wave  surface  of  the  extraordinary  ray 
is  a  (revolution)  ellipsoid.  To  be  able  to  apply  SNELLIUS' 
law  of  refraction  for  this  extraordinary  ray,  the  term 
"wavelength  standard  velocity"  is  introduced  instead  of 
"ray  velocity",  which  is  no  longer  equal  to  the  amount 
of  the  radius  of  the  (revolution)  ellipsoid  but  equal  to 
the  amount  of  the  radius  of  the  pertinent  (revolution) 
pedaloid. 
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In  /IS/  formula  (2.7)  Is  quoted  and  the  surface  is 
termed  * (revolution)  ovaloid" ,  the  pedal  locus  curve 
of  the  ellipse  is  termed  "oval",  since  in  the  above 
mentioned  case  of  double  refraction,  the  axes  of  the 
ellipse  etc.  differ  only  slightly.  The  shape  of  the 
pedal  locus  curve  of  the  ellipse  deviates  only  slightly 
from  that  of  an  ellipse,  that  is  why  this  term  ovaloid 
is  used. 

Special  case  2: 

sphere  -  if  all  three  axes  are  equal,  i.e.,  a  =  b  =  c, 
both  the  ellipsoid  and  the  tri axial  pedaloid  become 
spheres. 

If  the  passage  to  the  limit  c  -*•  0  is  affected,  the 
ellipsoid  degenerated  to  its  main  sectional  ellipse. 


2.6  The  (true)  pedaloid 

When  keeping  to  this  systematic  mathematical  classifica¬ 
tion  the  receiving  characteristic  in  the  elliptically 
polarized  field  can  be  interpreted  as  being  the  degenerate 
case  of  a  triaxial  pedaloid  by  setting  the  third  axis 
c  =  0.  In  contrast  to  the  ellipsoid,  the  three-dimen¬ 
sional  character  of  the  formation  thereby  developing  is 
preserved  in  this  case. 

The  parameter-free  representation  can  be  obtained  quite 
similarly  to  that  of  the  ellipsoid.  The  arbitrary 
radius  rpe  results  from  the  radius  of  the  pedal  locus 
curve  rp  (cf.  Eq.  (2.3))  by  multiplying  by  the  cosine 
of  the  angle  of  elevation  3 


rpe  “  rp  cos&  *  /a2cos2a  +  b2sin2a  *  cos  3  (2.8) 


By  substituting  the  angular  functions  as  mentioned 
before  the  equation  of  the  (true)  pedaloid 


a2x2  +  b2y2  =  (x2  +  y2  +  z2)2  (2.9) 

is  obtained  which  is  most  important  for  measurements  in 
the  elliptically  polarized  field.  The  families  of 
pedaloids  are  listed  in  Fig.  2.4. 

Special  case: 

Circular  ring  (torus)  with  the  inner  radius  being  zero; 
occurs  with  circular  polarization  of  the  field  (a  =  b) . 
Eq.  (2.9)  then  reads 

a2(x2  +  y2)  =  (x2  +  y2  +  z2)2  (2.10) 

The  arbitrary  radius  rt0r 

rtor  *  a  *  cos  8  (2.11) 


Degenerate  case: 

Double  sphere  -  representing  the  three-dimensional  re¬ 
ceiving  characteristic  in  the  linearly  polarized  field. 
On  account  of  b  =  0,  its  equation  reads 

ax  =  x2  +  y?  +  z2  (2.12) 

The  arbitrary  radius 

rds  =  a  cos  a  '  cos  8  (2.13) 


2.7  Intersection  of  the  pedaloid  with  a  plane 

Owing  to  the  measuring  technique  it  is  important  to  know 
the  sections  of  the  pedaloid  involved.  To  find  them,  a 
new  coordinate  system  x* ,  y' ,  z'  is  introduced  resulting 
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from  the  original  one  by  successive  rotation  by  the 
angles  u,  v,  w  (Eulerian  angles) .  The  origin  remains 
constant.  Setting  z'  =  O,  we  obtair.  the  cross-sectional 
area  lying  in  the  plane  <x',  y * ) . 

The  following  equation  describes  this  rotation 


r  =  (x,  y ,  z) 
r'  =  (x',y',z') 


(2.14) 


r  and  r'  are  position  vectors.  The  elements  of  the  matrix 
a^  are  the  direction  cosines  and  read 


aU 

" 

COS  V  • 

cos  w 

a  12 

= 

-cos  u 

•  sin  w 

a  13 

= 

sin  v 

a2l 

= 

cos  u  • 

sin  w  + 

sinu  • 

sin  v  • 

sin 

w 

a  22 

= 

cos  u  • 

cos  w  - 

sin  u  • 

sin  v  • 

sin 

w 

a2  3 

= 

-sin  u 

.  cos  V 

a  3 1 

= 

sin  u  ' 

sin  w  - 

cos  u  • 

sinv  • 

cos 

w 

a  32 

= 

sin  u  • 

cos  w  + 

cos  u  • 

sin  v  • 

sin 

w 

a  33 

= 

cos  u  • 

cos  V 

(2.15) 


In  Eq.  (2.9)  x,  y,  z  are  substituted  by  x' ,  y' ,  z'.  On 
account  of  the  above  relation  the  pedaloid  in  the  dashed 
coordinates  is  obtained: 


a2(a11x'  +  a12y'  +  a13z')2  +  b2(a2i*'  +  a22y'  +  a23z')2  = 
“  [(aux'  +  a12y 1  +  a13z')2  +  (a2ix*  +  a22y'  +  a23z')2  + 

<a3i  +  a32y 1  +  a33  z  ' )  2]  2 


(2.16) 
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Since  only  the  intersection  is  of  interest  here,  z'  =  0. 

By  squaring  the  brackets  in  Eq.  (2.16)  we  obtain 

a2  (ai  lX' 2  +  2a11a12x,y  +  a^y'2)  +  b^a^x'2  + 

+  2a21a22x'y'  +  a|2 y'2)  =  |  (a2 *  +  a21  +  a^Jx'2  + 

+  (a2r  +  af  2  +  a22)y ‘ 2  +  2 (aj  ja12  +  a2  xa2z  + a3ia32)x‘y ' 1 

(2.17) 

If  the  left  side  is  rearranged  and  if  we  further  consider 
on  the  right  side  that  the  sum  of  the  squares  of  the 
direction  cosines  is  equal  to  1  and  the  sum  of  the  com¬ 
posite  products  is  equal  to  zero,  we  obtain 

o  7  7  7  7  7  7  7  7  7 

(aaf,  +  b  a21)x'  +  (a  a12  +  b  a22)y'  + 

+  2(a2aMa12  +  b2a? ! a2 2 )x'y ’  =  U’2  +  y,?)2  (2.18) 

This  equation  again  represents  a  pedal  locus  curve  of  an 
ellipse  in  relation  to  the  coordinate  system  x',  y'.  On 
account  of  the  term  with  the  composite  elements,  major 
and  minor  axes  of  the  pedal  locus  curve  are  not  lying  in 
the  direction  of  the  coordinate  axes.  The  sections  on 
the  x'  ar.<2  y*  axes  are 

A'  =  / a2 a2  L  +  bza2]  and  B'  =  /a2a52  +  bza|2  (2.19) 

Proceeding  from  the  general  equation  of  sections  (2.18) 
the  various  special  cases  can  be  derived. 

Special  case  1;  (plane  vertical  to  main  pedal  locus  curve) 

a)  plane  vertical  to  the  original  pedal  locus  plane 
through  the  major  axis  a,  u  *  90°,  v  *  0° 
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Note:  Rotation  by  the  two  angles  u  and  v  is  sufficient 

to  reach  any  desired  direction  in  space.  The  matrix 
of  rotation  is  simplified  to  a  2  x  2  matrix:  the 
angular  functions  for  w  =  0  must  be  equal  to  1  or  0. 


The  matrix  of  rotation  therefore  reads 


aik 


1  O 
0  0 


Eq.  (2.18)  is  reduced  to 


(2.20) 


ax'  =  x'2  +  y'2  or 


y'2  =  ax'  -  x'2  (2.21) 


which  is  well  known  to  be  the  vertex  equation  of  the 
circle  with  the  radius  a/2,  where  x'  is  pointing  into 
the  direction  of  x  and  y*  into  that  of  z. 


b)  plane  vertical  through  minor  axis  b,  u  =  0°,  v  =  90° 


matrix 


aik 


0  0 
0  1 


(2.22) 


In  this  case, 

by'  =  x'2  +  y,z  or  x'2  =  by*  -  y*2  (2.23) 

is  the  vertex  equation  of  the  circle  with  the  radius  b/2, 
where  x'  is  pointing  into  the  direction  of  z  and  y* 
into  that  of  y. 


c)  arbitrary  plane  vertical  to  main  pedal  locus  curve, 
u  =  90°,  0  <  v  <  90° 


matrix 


aik 


cos  V 

0 

sin  v 

0 

(2.24) 


yields 


r  •  2 


X 


•  2 


x*  /  a2cos2v  +  b2sin2v 


(2.25) 
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The  radical  represents  the  radius  of  the  intersected 
pedal  locus  curve  and  is  equal  to  the  diameter  of  the 
vertex  circle  represented  by  Eq.  (2.25). 

Result  from  1  a)  -  c) :  All  surfaces  of  intersection 
vertical  to  the  main  pedal  locus  curve  of  the 
pedaloid  are  vertex  circles  in  relation  to  the 
center  of  the  pedaloid. 


Special  case  2: 

Intersections  through  the  minor  axis  b,  u  =  0°,  0<v<  90° 

| cos  v  0 
0  1 

We  obtain  the  pedal  locus  equations 


matrix 


aik 


(2.26) 


aV  x'2  +  by'  =  (x*2  +  y’2)2  (2.27) 

u 

The  condition  an  =  cos  v  =  b/a  yields  a  circle  with 
the  radius  b  as  surface  of  intersection. 


Proof: 


a2  x' 2  +  b"y 1 2  * 

az 

(x'2 

+  y'2)2 

b2 (x'2  +  y,? )  * 

(x'2 

+  y ,2 ) 2 

x 1 ~  +  y ' 2  *  b 


K  o  ? 

If  cua  v  ■  —  ,  b  >  a2 an  and  forms  the  major  axis  of 

cl 

the  curve  of  intersection  which  is  again  a  pedal  locus 


curve. 
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Special  case  3: 

Intersections  through  the  major  axis  a,  0  <  u  <  90° ,  v  =  0° 


matrix 


ik 


1  0 
0  cos  u 


(2.28) 


We  obtain  the  pedal  locus  equation 


a2x*2  +  b2a|2y'2  -  (x'2  +  y'2)2 


(2.29) 


Summarizing  the  following  characteristic  properties  of 
the  pedaloid  can  be  given: 

1.  The  main  intersection  (containing  the  axes  a  and  b 
of  the  ellipse)  is  the  pedal  locus  curve  with  the 
largest  area. 

2.  The  radius  vector  vertical  through  the  center  0  of 
the  main  intersection  is  equal  to  zero. 

3.  All  intersections  vertical  to  the  main  pedal  locus 
surface  have  the  shape  of  a  figure  eight.  The  largest 
area  of  this  kind  occurs  at  the  intersection  con¬ 
taining  the  maj or  axis  a,  the  smallest  area  arises 

at  the  intersection  containing  the  minor  axis  b. 

4.  That  intersection  area  that  contains  the  minor  axis 
and  is  inclined  to  the  main  pedal  locus  curve  by 
the  angle  v,  so  that 

cos  v  =  a/b 

is  valid,  is  a  circle  with  the  radius  rpe  =  b. 


2.8  Calculation  of  the  squint  at  polarization  ellipses 
lying  arbitrarily  in  space 

If  the  polarization  ellipse  is  not  lying  in  the  horizontal 
plane  of  measurement  and  if  the  direction  finder  for  the 
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detection  of  the  maximum  : A  minimum  t£  field  strength 
is  rotn:.d  only  in  this  plane,  a  specific  error  of 
measurer-.  :.t  arises  which  can  be  determined  by  the  fol¬ 
lowing  cm  rome  value  problem. 

If  the  coordinate  system  in  which  the  major  axis  of  the 
pedalcid  is  pointing  in  the  direction  of  the  x-axis  and 
the  minor  axis  in  y-axis  direction,  is  revolved  on  the 
u-axis  by  the  angle  u  and  consequently  on  the  new 
y I -axis  by  the  angle  v,  any  arbitrary  intersection  area 
can  be  obtained  as  x'y '-plane.  Let  this  intersection  arc.-, 
of  the  podaloid  be  the  pedal  locus  curve  measured  hori¬ 
zontally;  the  question  aiises,  what  si2e  the  angle  v  iaust 
assume  ti>  make  the  newly  developed  axes  x'  and  y'  point 
in  the  directions  of  a*  and  b',  respectively;  a'  and  b* 
are  the  major  and  minor  axes  of  the  pedal  locus  curve  in 
the  intersection.  This  angle  ■  is  at  the  same  time  that 
squint  that  is  made  with  the  above  described  direction 
finding . 

Proceeding  from  the  equation  of  the  pedaloid  in  the  dan  ..  . 
syster ,  the  coef  icient  of  x'  living  the  length  of  th>> 
section  A  *  on  the  x'-axts  reads  as  fellows: 

A  *  ’  &  aj  5  +  b  a  . . )  - 

-•  a  cos  v  «  cos  w  +  b;  (cos?u  ♦  sin^w  + 

•  •> 

'•  sin  u  sin  v  cr. s' w  +sinu  •  cos  u*  sin  v  •  sin  7 

(2. 300 

The  derivation  with  respect  to  w  is  formed  and  set  zero 
resulting 

sin  2u  *  sin  v 

-  tan  2  =  — - 

—  cos'v  +  sin2u  •  sin2v  -cos;u 
b7 


sin  2w 
cos  2w 


< 2  .  V  ! : 


BEARING  ERROR  W 


FIG.  2.6 
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For  a  check  the  angle  obtained  from  the  above  equation 
can  be  inserted  into  the  coefficient  of  the  "term  of 
rotation"  of  Eq.  (2.18)  which  then  has  to  vanish.  The 
easier  way  is  to  set  the  coefficient  C'  of  the  term  of 
rotation  equal  to  zero  and  calculate  w  therefrom  and 
compare  it  with  the  above  result  proving  it. 

The  angle  w  was  consequently  computed  for  the  following 
parameters : 

0  (=  a  :  b)  =  1;  1.5;  2;  2.5;  3;  3.5;  4;  4.5;  5 

“j-  =  1,  2,  3,  4,  5,  6,  7,  8,  9,  10 ,  20,  30,  ...,  80° 

The  results  are  represented  graphically  in  Figs.  2.5  >2.8 
for  different  ratios  of  axis  a  :  b  of  the  polarization 
ellipse.  Fig.  2.5  is  valid  for  a  circularly  polarized 
field  (a  =  b) .  From  the  other  diagrams  it  is  to  be  seen 
that  only  for  great  angles  u  and  v  the  squint  is  essen¬ 
tial,  e.g.  for  u  =  v  =  30°,  w  =  12°  (a  :  b  =  3  :  2) , 
cf.  Fig.  2.6.  For  ellipses  with  a  smaller  ratio  of  axes 
this  squint  becomes  less  important.  Fig.  2.8  shows  that 
a  distortion  of  the  ellipse  u  =  v  =  30°  yields  the 
negligible  angle  w  =  4°. 


2.9  Conclusion-: 

With  the  mathematical  representation  of  the  receiving  cha¬ 
racteristic  in  an  elliptically  polarized  electromagnetic 
field  (=  pedaloid)  it  is  possible  to  calculate  every  arbitrary 
section  plane  of  this  pedaloid. inis  fact  involved  many  ad¬ 
vantages  for  field  strength  measurements.  One  of  the  most 
important  things  is  that  we  are  able  to  calculate  the  squint 
that  arises  when  the  direction  finder  is  rotated  only  in 
the  horizontal  plane  when  measuring  the  maximum,  the  pola¬ 
rization  ellipse,  however,  lies  arbitrarily  in  space. 


CHAPTER  3 


A  NEW  METHOD  TO  DETERMINE  THE  SPATIAL  POSITION 
OF  THE  MAGNETIC  ROTARY  FIELD  ELLIPSE 


3.1  Introduction 

For  many  experiments  concerning  VLF-reception  it  is 
necessary  to  know  the  exact  position  cf  the  rotary  field 
ellipse.  N.  NESSLER  described  a  method  to  determine  the 
pedaloid  in  his  thesis  /10/.  This  method,  however,  showed 
to  be  time-wasting  and  in  some  cases  unexact;  the  error 
made  by  this  method  is  calculated  by  R.  Klein  /16,  17/. 
The  problem  was  to  find  a  new  method  using  all  charac¬ 
teristics  of  the  pedaloid  for  a  better  and  non-sub jective 
calculation.  All  the  quantities  of  the  rotary  field 
ellipse  had  to  be  measured  by  means  of  an  antenna  system 
constant  in  space  instead  of  using  the  direction  finding 
antenna  of  former  experiments. 


3 . 2  The  pedaloid 

The  measurement  of  a  rotary  field  ellipse  always  leads 
to  "pedal  locus  curves"  when  using  ferrite-rod  direction 
finding  antennas. 

Let  us  now  briefly  consider  this  effect  mathematically: 

a  ...  maximum  voltage  -  major  axis 
b  ...  minimum  voltage  -  minor  axis 
u  ...  angular  frequency 

4>  ...  angle  of  inclination  between  major  axis 
and  antenna 

...  angle  between  the  normal  of  the  ellipse 
plane  and  the  antenna 


0 


The  voltage  Induced  into  an  antenna  in  dependence  on  u>, 
4 ,  and  0  is  given  by 

U  =  (a*cos<f> -cosut  +  b*sin$ ’sinwt) sin9  (3.1) 


with  a  and  b  containing  all  conversion  factors  between 
magnetic  field  strength  and  voltage. 

The  temporal  maximum  (minimum)  of  the  voltage  is  reached 
when 


dU 

dt 


=  0  =  u  (-a*cosiji  *sinut0 +b*sin^ ‘CoscDto)  sinB 


t-tQ 


(3.2) 


resulting 


& 


tan$  =  tanutQ 


(3.3) 


Which  further  on  will  be  referred  to  as  "pedaloid- 
condition" . 


It  can  be  seen  that  the  phase  angle  relative  to  the 
major  axis  is  not  identical  with  the  spatial  angle, 
except  for  the  case  a  =  b.  And  the  antenna  does  not 
reproduce  the  ellipse,  but  a  pedal  locus  curve  or  a 
pedaloid . 

3 . 3  Theory  of  the  new  method 

The  antenna  combination  consists  of  three  electrically 
identical  ferrite-rod  antennas  perpendicular  to  one 
another  thereby  forming  an  orthogonal  trihedral.  The 
three  axes  of  the  pedaloid  are  represented  by  the  vectors 
eaf  ®b»  and  n,  with  the  subscripts  a  and  b  standing  for 
the  major  and  minor  axes  of  the  ellipse,  respectively, 
and  n  being  the  normal  vector  of  the  rotary  field  ellipse. 
The  arbitrary  position  of  the  rotary  field  ellipse  in 
respect  to  an  other  Cartesian  system  is  given  by  the 
three  Eulerian  angles  $,  1 1>,  0.  Besides,  two  more  unknown 
quantities  are  requested  for  a  complete  description  of 
the  ellipse  In  space,  namely  the  two  main  voltages  of 
the  ellipse.  If  we  succeed  in  measuring  these  five 
quantities  position  and  size  of  the  ellipse  arbitrary 
in  space  are  definitely  determined. 

3.4  Calculation 

Fig.  2  shows  the  normal  of  the  ellipse  in  reference  to 
the  &:tenna  system. 

Let  the  directions  of  the  antennas  be  i,  j,  and  k, 
the  following  conditions  can  be  calculated  (antenna 
directions  expressed  by  direction  cosines  in  respect 
to  ea,  e^,  and  n) : 
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Pig.  3.2 


H-t 

II 

cos^j • sin6 j 
sin41 -sine , 

il 

to 

•  «* 

cos  if  2  •sin02' 
sinf 2  *8in82 

II 

+M 

cosf 3 • sin6 3 ' 
sin<t>3  »sin83 

1 

cos9  i 

1  1 

COS02 

1  1 

cose  3 

(3.4) 

■>  *►  •¥ 

As  i,  j,  and  k  form  an  orthogonal  trihedral  the  following 
conditions  are  valid,  too. 

t  t 

i  *  j  ~  k 

j  *  k  =  i  (3.5) 

k  *  i  —  j 

Together  with  the  pedaloid  condition  (3.3)  we  get  from 
(3.1)  the  voltages  in  terms  of  <J>j_  and  6^: 

=  (a2*cos2i}>i  +  b2  *sin2<{ii) •  sine^  (3.6) 

and  with  the  abbreviation  u>t0  =  il  we  obtain 


II 


A 
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Ui  =  a*b*  (a2*sin2fii  +  b2  •cos2fli)-1^2  •  sindj^  (3.7) 

and  thereby  the  voltages  expressed  by  the  phase  angles  0^. 
These  three  voltages  form  a  vector  with  a  constant  amount , 

e.g. 

3 

J  Uj2  =  a2  (cos241  ‘Sin^!  +  cos242  *sin202  +  cos243*sino3)  + 
i=l 

b2  (sin2^  *sin2 6j  +  sin242 *sin2e2  +  sin243 *sin03) 

(3.8) 

From  Eq.  (3.4)  we  know  that  cos^-sinBjL  and  sin4i*sir.0i 
are  the  direction  cosines  of  the  pedaloid  axes  ea  and  e^ 
in  reference  to  the  antenna  system.  As  the  sum  of  the 
squares  of  direction  cosines  always  equals  1  we  may  write 
formula  (3.8)  shortly  as 

3 

l  U.2  =  a2  +  b2  (3.9) 

i=l  1 

On  the  other  hand,  the  condition  of  orthogonality  de¬ 
livers  the  following  formula: 

k  *  i  x  j  = 

*2^3  ”  i3^2  | 

i33 i  ”  ii33  j 

_  (cos^j  *sin42  **  cos$2  *sin<t>  j )  sine  l  •  sine ?  j 

(3.10) 

The  main  result  of  this  formula  is: 

cose 3  =  -sin(<)»i  -  <?2)sin6i  *sin02 
cose 2  =  -sin($2  “  ^iJsinej *sin03 

COS0J  =  -sin(<|>3  -  4»3) sine2  •  sines 
With  the  pedaloid  condition  (3.3)  we  derive 

sinUj  -  4> 2 )  =  simfrj .cos$2  -  sin<fr2  *cos<J>  = 


r 

^COS0  3 


(3.11) 
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We  then  obtain  with  Eq.  (3.7) 

sin(Q,  -  fl,) 

sint*!  -  *2)  =  -  l^Uj,  (3.13) 

ab'SinOj  *sine2 

Consequently,  Eq.  (3.11)  with  the  abbreviations 
fij  -  fl2  =  Xj ,  fl2  -  J)3  «  X2 ,  &3  -  ni  =  X3  (phase  dif¬ 
ferences  between  two  antennas)  reads  as  follows: 


COS0  ! 

sinX2 

ab 

U2U3 

ft 

sinX3 

[T  11. 

^US  D  £ 

ab 

U3U1 

sinXx 

IV  ft 

COS  8  ^ 

ab 

UiU2 

Since  these  cosines  cose^  are  the  direction  cosines  of 
the  normal  of  the  ellipse  in  reference  to  the  antenna 
system 

l  cos2^i  =  1  (3.15) 

i=l 

is  valid,  or,  for  further  calculations 

a2b2  =  sin2X1U^U2  +  si.n2X2U2U2  +  sin2X3U3uf  (3.16) 

Formulas  (3.9)  and  (3.16)  together  allow  the  evaluation 
of  a  and  b. 
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P2  =  a2  +  b2  =  l  Ui2 
iii  1 


(3.17) 


Q2  =  a2  •  b2  =  f  sin2X.O,2 •  D, 
i=l  11  J 


(3.18) 


and  there from 


a2  .  *L  w  _El  .  q2 

2  4 


(3.19) 


(3.20) 


So  the  two  main  voltages  are  known  in  texms  of  the 
measured  quantities! 


3.5  Calculation  of  the  Eulerian  angles 

The  three  Eulerian  angles  are  0,  and  $. 

e  ...  angle  of  nutation  (between  ellipse  normal  n 

and  the  direction  k) 

4/  ...  angle  of  precession  (between  axis  i  and  inters*?  ion 

line  x  of  the  ea-ejj -plane 
the  i-j-plane) 

4>  ...  angle  of  torsion  (between  the  intersection  line  x 

and  the  axis  e^) 


0:  It  is  evident  (Fig.  3.2)  that 


=  e3 


(3.21) 


ifr:  First  we  have  to  determine  the  intersection  line  x. 
So  we  fc-rm  the  cross  product  of  the  vectors  ic  and  n. 


sin<t>3 

k  x  u  =  -cos4>  3  =  x 


(3.22) 
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Mow  the  req’V'ated  angles  are  easily  to  be  calculated: 

(X  •  i)  =  cos*  =  cos* 1  * sin6  j  *  sin*3  -  sin* j  •  sine  j  *cos*3 

(3.23) 

As  there  are  certain  conditions  valid  between  the  direction 
cosines  in  the  matrix  of  torsion  we  can  write  Eq.  (3.23) 
as  follows: 

->  +  cose. 

(X  •  i)  *  cos*  - -  (3.24) 

sine  3 

Together  with  (3.15) 

cose, 

sin*  =  —  —  -  (3.25) 

sine3 

follows . 


*:  The  angle  *  follows  from  the  scalar  product  of  the 

-*■  ,  4 

vectors  x  and  ea. 

x  •  ea  =  cos*  =  sin*3  (3.26) 

sin*  follows  from  the  cross  product  of  these  two  vectors: 


lx  x  e. 


=  sin*  =  cos*. 


(3.27) 


We  now  have  to  write  the  Eulerian  angles  in  terms  of  the 
voltages  Uj ,  U,,  ,  U3  and  the  phase  angles  Xx  ,  X2  ,  X3  t 


6: 


sinX 


cose  5  cos ft  =  - 


l 


ab 


UiU2 


(3.28) 


As  ab  is  evaluated  in  Eq.  (3.16),  cose  is  known. 


±i 


cos*  = 


COS02 

Sin03 


(3.29) 
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sin  3 

cose,  »  -  — -— =■  0,0, 
*■  ab  3  1 


(3.30) 


sinij  »  A  -  cosze3  *  ^  •  /sin2X20|of  +  8in2X3oiOi2  (3.31) 


COS* 


f- 


sin2X2022 


sin2X3U2 


+  1 


(3.32) 


sin*  * 


COS  01 


CO80J 

sine  3 

sinX1U1U2 

ab 


(3.33) 


(3.34) 


sin* 


-1 


/  .  sin2X30!2 

r1  +  - — T 

sin4X2U22 


(3.35) 


So  the  angle  *  Is  known,  too. 


•Is  From  Eq.  (3.6)  we  get: 


cos*  *  sin*-j 

(3.36) 

O32  -  (a^cos2^^  +  b2sin2*3)sin203 

(3.37) 

^  ... 

?  sln?e3 

sin7« ,  *=■  - -  - — 

3  b2  -  a2 

(3.38) 

sin*  b  cos*3 

(3.39) 

U32 

-  b2 

sin2e, 

cos  •}■'  •>  =  - 5-2 - 5- 

a2  -  b2 

(3.40) 
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As  a,  b,  and  sin6g  are  known  from  Eqs.  (3.19),  (3.20), 
(3.14),  all  Eulerian  angles  are  known  In  terms  of  the 
measured  quantities. 


3.6  Conclusion 


It  is  possible  to  determine  the  exact  spatial  position 
of  the  rotary  field  ellipse  by  measuring  three  voltages 
and  three  phase  angles .  The  results  can  be  computed. 

This  new  method  has  the  advantage  of  being  quick  and  free 
of  errors  of  adjustment  (measurement  possible  at  any 
trihedral  positionl)  and  can  be  done  automatically  with¬ 
out  any  further  adjustment  of  the  antenna  position. 


CHAPTER  4 

EXTERNALLY  KEYED  PHASE-LOCK  DEVICE 


4 . 1  Equipment 

The  externally  keyed  (by  VLF-wave  propagation)  phase- 
lock  device  described  in  /ll/  measures  the  phase  dif¬ 
ference  between  a  VLF  wave  with  the  frequency  of 
2.667  kHz  in  the  near  field  and  the  GBR  signal  (16.0  kHz) 
in  the  far  field.  The  measuring  principle  shall  be  brief¬ 
ly  repeated  (cf.  Fig.  4.1,  4.2,  4.3). 

The  GBR  signal  (A2  or  F2  keyed)  is  received  by  ferrite 
rod  antennas  and  used  for  phase-stable  triggering  of 
an  extremely  stable  balance  oscillator.  The  lock-in 
range  is  only  20  Hz  so  that  only  the  16,000  Hz-signal 
is  used  for  phase  triggering  during  the  FM-keyed  time 
(1C, 000  and  15,950  Hz),  thus  avoiding  influences  of 
keying  effects  and  giving  a  stable,  extremely  exact  time 
normal  at  any  arbitrary  point  of  reception. 

This  type  of  time  normal  is  used  to  trigger  both  the 
transmitter  and  the  receiver  of  a  VLF-signal  transmis¬ 
sion  chain.  Before  triggering  the  frequency  must  be 
changed  in  order  not  to  affect  the  original  "clock¬ 
triggering"  signal.  To  get  the  required  high  accuracy 
only  a  digital  frequency  change  is  allowed.  In  the 
described  equipment  a  divide-by-6-mode  is  used  so  that 
the  frequency  of  2.667  Hz  results. 

In  the  transmitting  equipment  this  signal  is  amplified 
and  sent  by  a  frame  antenna. 


The  receiving  equipment  uses  two  ferrite  rod  antennas 
for  the  reception  of  the  VLF  signal  and  of  the  GBR  syn¬ 
chronous  signal ,  respectively.  The  two  signals  (VLF 
signal  and  the  divided  16  kHz  signal)  are  compared  in 
the  phase  detector.  The  result  is  a  DC-signal  changing 
with  phase  difference  which  may  act  as  phase  propor¬ 
tional  voltage  on  the  one  hand  or  as  rectified  signal 
of  a  receiver  with  extremely  small  bandwidth  (depending 
on  integration  time  constant;  here  0.2  Hz!)  on  the 
other  hand. 

The  measuring  practice  showed  the  necessity  of  some 
changes  and  supplements  in  circuit  design,  e.g.: 

1)  a  special  Schmitt -trigger  to  make  the  triggering 
correspond  with  the  zero-axis  crossing  of  the  GBR 
signal,  thus  ensuring  that  the  phase  relation  is  inde¬ 
pendent  on  the  amplitude  (further  details  see  Fig. 

4.1  or  /ll/) i 

2)  a  measuring  amplifier  for  measuring  the  incoming 

VLF  signal  combined  with  an  attenuator  providing  a  uniqu 
signal  amplitude  at  the  phase  detector  independent  on 
the  input  voltage  (cf.  Fig.  4.2  and  4,3); 

3)  a  main  operation  mode  switch  with  the  positions 

"off  /  battery  control  /  GBR  amplitude  /  VLF  amplitude  / 
phase" . 

A  number  of  test  jacks  were  added  at  the  front  panels 
of  both  the  transmitting  and  receiving  parts  for  opera¬ 
tion  control  and  the  possibility  of  connecting  external 
devices  (e.g.  line  recorder) . 

The  complete  electrical  and  mechanical  design  is  shown 
in  Figs.  C,  D,  E,  F. 


4B 


4.2  iiea.-arements 


The  a.l-,j.t"u.nts  were  scheduled  for  the  first  quarter 
of  it  could  net  be  done  in  time  because  GBR 

stopx-c  ;  . raasmitting  in  January  for  a  longer  period 
of  t  i  . 

f*.-: .•’»n  had  starts 3  again ,  the  phase  measurements 
veto  ....... '..tod  in  St.  C-^rtraudi  mine  /  Tyrol.  The 

ill*; ^  s.jriet  proveu  to  be  correct  because  it  was 

r c:-.e  • .  fwice. 

The  amplifier  (1  kW  maximum  output)  was  used 

as  ■  ctor  in  cennecti-.n  with  the  big  loop  (40x40  ir.. 


TO  i'  .. r. }  which  was  tu;.^d  to  resonance  by  capacitors. 
Ait!:...*  ;  -'sore  is  a  stepper  circuit  for  the  VLF  fre- 
-i'i .  \  .  uircnt  in  the  antenna  (15  A)  produces  such 

j  i  v;  field  in  -he  very  near  measuring  points, 

c.la:  ferrite  red  antenrsc*  for  the  reception  of  the 

**oi:  t  •  *  .w  CBP  si  gnal  hod  to  be  oriented  very  exactly 
in  i.  field  strength  direction  of  the  VLF  signal 

in  c: _c-  prevent  an  overload  of  the  input. 

The  :  e  ing  points  were  situated  along  the  "West- 
seals g ’  - -:-?llery  which  is  nearly  parallel  to  the  Inn 
valle.  '  Le  4.1  gives  the  measuring  results. 
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Table  4.1 


d 

|m| 

UVLF 

|nV| 

Phase 

l°l 

100 

27 

20 

150 

21 

15 

200 

14 

15 

250 

7.4 

15 

300 

5.0 

20 

350 

2.5 

20 

400 

1.9 

20 

450 

1.5 

25 

500 

1.1 

30 

550 

0.9 

20 

600 

0.77 

18 

650 

0.63 

22 

700 

0.54 

15 

750 

0.47 

0 

distance  between  transmitter  and  receiver 
voltage  proportional  to  received  VLF  signal 
phase  difference  between  synchronization 
signal  GBR  and  transmitted  VLF  signal  reduced 
modulo  60°  (because  of  the  6  :  1  reduction  of 
the  GBR  frequency) 

Result:  Primarily  it  was  assumed  that  the  GBR  incoming 
direction  is  vertical  to  the  line  transmitter  -  re¬ 
ceiver  so  that  there  is  no  phase  difference  between 
these  two  points.  Therefore  thi  difference  was  ex¬ 
pected  to  rise  with  increasing  distance. 


d . 

Uvlf  •  • 

Phase  . 
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4 . 1  Co n  c 1 us ion 


Th?o  phase  measuring  device  measures  the  phase  dif- 
fexencc  between  the  synchronization  signal  and  the 
VLF  signal,  and  has  the  disadvantage  that  for  measuring 
the  phase  o£  one  channel  the  phase  shift  of  the  other 
-r. annex  «.ust  be  known-.  If-transmitter  and  receiver  are 
exchanged  the  phase  shi fts  between  VLF  and  GBR  signal 
can  be  added  and  subtracted,  respectively.  Thus,  it  is 
possible  ro  eliminate  the  phase  shift  of  the  GBR  signal. 


CHAPTER  5 


A  DIGITAL  PHASE  MEASURING  DEVICE 


5 . 1  Introduction 

To  be  able  to  measure  phases  at  the  propagation  of 
electromagnetic  waves  in  the  VLF  range  a  device  was 
developed  by  Dr.  N.  NESSLER  /20/  showing  the  possibi¬ 
lity  of  phase  measurements  within  this  frequency  range 
for  the  first  For  the  synchronization  between 

transmitter  ar.u  receiver  a  further  external  transmitter 
was  used.  This  existing  equipment  allows  phase  measure¬ 
ments  only  at  a  certain  frequency,  and  in  some  special 
cases  indication  is  amplitude -dependent.  A  new  phase 
measuring  device  comparing  two  incoming  signals  (measuring 
signal  and  reference  signal)  completely  independent  on 
frequency  and  amplitude  is  now  being  built. 

The  requirements  which  should  be  fulfilled  by  this  de¬ 
vice  are  rather  high:  current  consumption  must  be  very 
small  to  enable  battery  operation;  the  device  should 
have  little  weight  and  volume  and  be  insensitive  to  heavy 
pushes  and  high  atmospheric  humidity  to  be  able  to  be 
operated  in  the  mine;  production  cost  must  be  kept  low 
(since  the  main  part  of  the  device  is  intended  to  be 
used  as  well  for  a  special  measuring  equipment  of  the 
rotary  field  ellipse) , 

Thorough  information  showed  that  no  commercial  device 
fulfills  all  these  requirements ,  that  is  why  it  was  con¬ 
sidered  most  efficient  to  build  it  in  our  own  laboratory. 

The  device  is  to  indicate  the  phase  digitally  independent¬ 
ly  on  amplitude  and  frequency  within  the  frequency  range 
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of  100  Hz  -  100  kHz  over  a  liquid  crystal  display 
counting  from  -180°  to  +180°. 


5.2  Basic  construction 

Fig.  5.1  shows  the  block  diagram  of  the  device  without 
a  representation  of  receiving  antennas  and  preamplifiers. 
The  two  receiving  signals  the  mutual  phase  position  of 
which  is  to  be  measured  are  connected  to  the  inputs  A 
and  B.  The  two  Schmitt  trigger  units  transform  the  sine 
signal  at  the  input  to  a  square-wave  signal  with  the 
pulse  duty  factor  1  :  2  (pulse  duration  :  pulse  period) . 
The  two  square-wave  signals  are  compared  as  to  their 
coincidence  in  the  discriminator.  The  pulse  duration 
of  the  output  signal  at  the  discriminator  is  a  direct 
measure  for  the  phase  shift.  In  the  logic  unity  this 
phase-proportional  signal  cuts  in  a  counting  oscillator 
(10  MHz)  the  number  of  pulses  of  which  is  counted.  This 
number  within  a  definite  time  is  a  direct  measure  for  the 
phase  angle  independent  on  input  amplitude  and  input  fre¬ 
quency  and  is  indicated  at  che  liquid  crystal  display. 


5.3  Description  of  circuit 

5.3.1  Input  Schmitt  trigger  with  null  hysteresis 

Fig.  5.2  shows  the  circuit  for  channel  A.  The  sine  signal 
is  connected  to  input  A  (anout  30  mVpp  -  5  Vpp) .  In  the 
upper  branch  an  impedance  transformer  follows  (T  1)  and 
connected  thereto  the  one  "half-Schmitt  trigger".  In  the 
bottom  branch  the  sine  signal  is  inverted  by  a  T  2  and 
the  gain  of  the  transistor  is  set  to  1  by  a  P  3.  The 
succeeding  T  3  works  as  an  emitter  follower  and  controls 
the  second  half-Schmitt  trigger.  Such  a  half-Schmitt 
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FIG.  5.3 


through.  This  interruption  is  necessary  when  the  input- 
signal  becomes  discontinuous  as  it  is  the  case  with 
keyed  transmitters.  Counter  2  can  reach  360000  at  most 
since  the  two  counters  only  work  when  A  =  L  (pulse  duty 
factor  of  A  is  1  :  2) . 

After  dividing  this  number  by  two  and  cancellation  of 
the  last  three  decimal  places  the  content  of  the  counter 
corresponds  to  the  phase  shift  measured  in  angular 
degrees.  The  end  of  one  counting  cycle  sets  the  mono¬ 
stable  multivibrator  I  to  a  logical  'L'  which  blocks  the 
10  MHz  clock  pulse  over  NAND  gate  II.  At  the  same  time 
the  bistable  latches  are  switched  through.  The  first  two 
decimal  places  of  counter  2  and  the  carry  over  of  counter 
2  are  directly  transferred  to  the  indicator.  The  falling 
slope  of  MMV  I  switches  a  second  monostable  multivibrator 
(MMV  II)  to  the  quasistable  position.  Consequently  all 
counters  and  the  flip-flop  for  the  carry-over  (NAND  gates 
III  and  IV)  are  zeroed;  a  new  counting  cycle  can  start. 

The  above  described  arrangement  allows  a  phase  indica¬ 
tion  independent  of  the  input  frequency  and  independent 
of  the  10  MHz  clock  frequency. 

5.3.4  The  display  unit 

A  2  1/2-place  liquid  crystal  display  (L.C.D.)  is  used 
which  has  the  advantage  of  an  extremely  low  power  con¬ 
sumption  (about  40yu.W/ segment,  i.e.  less  than  1/1000  of 
the  power  requirement  of  a  segment  of  an  L.E.D.  The  drive 
of  the  L.C.D.  is  more  laborious  than  triggering  the  L.E.D. 
since  the  L.C.D.  can  only  be  operated  by  alternating 
current  from  about  15  Veff  at  20  -  100  Hz.  A  special 
circuit  was  now  developed  producing  an  a-c  voltage  of 
60  Vpp  which  is  shown  in  Fig.  5.5. 
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The  a-c  voltage  is  pro-iuced  by  a  self-starting  multi¬ 
vibrator  /22/  oscillating  at  a  frequency  of  about  80  Hz. 
The  three  coil  windings  were  bridged  by  capacitors  (C2, 
C3,  C4)  to  avoid  dangerous  voltage  peaks  at  the  output. 
The  output  voltage  becomes  sinusoidal  and  nearly  in¬ 
dependent  on  the  load.  This  output  voltage  is  now  recti¬ 
fied  over  C5  and  Dl.  With  this  circuit  (Villard-doubling) 
the  a-c-componeni  is  maintained  but  is  increased  by 
1/2  Vpp  (pulsati  g  d-c  voltage  between  0  V  and  60  V). 

The  transistors  are  special  high-blocking  types  with 
°CE0  -  60  V.  Only  the  a-c-component  is  coupled  out  over 
C6-C21  and  laid  onto  the  segment  of  the  L.C.D. 

By  using  this  circuit  we  can  operate  the  display  with 
the  required  a-c  voltage  of  about  20  Veff  and  80  Hz, 
which  grants  a  high  lifetime  of  the  segments.  At  the 
development  of  this  circuit  it  was  taken  into  account 
to  keep  power  consumption  low.  The  current  input  for  the 
whole  display  unit  is  not  more  than  3  mA  at  5  V. 

5.3.5  Power  supply 

The  device  is  operated  by  2  batteries.  For  the  linear 
circuits  a  12  V  accumulator  is  used.  The  voltage  of  i  6V 
for  the  operational  amplifier  is  obtained  by  voltage  di¬ 
vision.  The  logic  is  operated  by  a  separate  battery 
(Fig.  5.6).  An  estimation  of  the  power  consumption  when 
using  normal  ICs  showed,  that  the  logic  would  take  about 
3300  mW.  Consequently,  it  is  advisable  to  use  low-power 
ICs  wherever  possible  to  reduce  power  consumption  to 
about  one  half.  In  spite  of  the  high  prime  costs  the  use 
of  low-power  ICs  is  justified,  since  it  will  then  be 
possible  to  take  smaller  and  cheaper  batteries  and  to 
reduce  weight  and  volume  of  the  device.  Since  the 
switching  time  of  the  *L* -series  is  longer  than  that 
of  the  normal  series ,  low-power  ICs  can  only  be  used 
where  the  limiting  frequency  is  below  1  MHz. 
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5.4  Error  estimation 

The  kind  of  logic  causes  2  errors  to  arise  at  counting 
the  Q  signal. 

Error  1 

At  counting  the  Q  signal  by  means  of  the  10  MHz  signal 
the  number  of  clock  pulses  has  an  inaccuracy  of  ±1  pulse 
per  impulse  train  (Fig.  5.8) .  During  one  counting  cycle 

ni  *  *  vi 

t  ...  duration  of  a  counting  cycle  (sec) 
v1  ..  measuring  frequency  (Hz) 

of  such  impulse  trains  occur.  The  content  of  counter  2 
has  a  variation  of  t  m  in  the  most  unfavorable  case. 

Translated  to  angular  degrees  the  indication  error 
amounts  to 

*  ±t*v1 »360/a 

At j  ..  error  1  (°) 

a  ....  maximum  content  of  counter  1 

(in  this  case  fixed  to  360000) 

With 

t  =  a/v2 

v2  ...  clock  frequency  (Hz) 
we  have  a  final  result  for  error  1 

A<j>j  =  ±^^2*360  (5.1) 

This  shows  that  the  error  depends  only  on  the  ratio  of 
input  frequency  to  clock  frequency.  An  improvement  can 
be  reached  by  increasing  the  clock  frequency.  In  the 
present  device  v2  =  107  =  10  MHz. 
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As  is  shown  in  Table  5.2,  error  1  is  relatively  large. 

A  control  at  the  ready-built  device,  however,  showed 
the  real  error  to  be  smaller  than  in  Table  5.2  by  one 
order  of  magnitude,  which  is  due  to  the  fact  that  the 
excess  and  the  lack  of  counted  pulses  neutralize  each 
other.  An  exact  calculation  by  means  of  calculus  of 
probability  confirms  this  measurement  result. 

Error  2 

Fig.  9  shows  the  formation  of  error  2;  we  can  distinguish 
three  cases: 

Case  a)  is  the  ideal  case.  In  the  impulse  interval  of 
the  impulse  train  the  MMV  rebounces  to  the  stable  state 
thereby  starting  the  counting  cycle  as  soon  as  A  =  L. 

Let  the  ratio  of  clock  frequency  to  measuring  frequency 
be  a  wno] 3-number  multiple  modulo  360,000  of  frequency  A, 
so  that  after  a  definite  number  of  complete  impulses  A 
360,000  clock  impulses  were  counted  in  counter  1.  This 
means  that  every  complete  impulse  train  of  the  duration 
tl  corresponds  to  a  complete  impulse  train  t .  Since  the 
ratio  tl  :  c  which  is  proportional  to  the  phase  is 
maintained,  no  error  of  the  second  type  arises  in  this 
case. 

Case  b)  Let  the  frequency  of  A  be  somewhat  smaller  than 
for  case  a,  so  that  the  360,000  impulses  of  counter  1 
do  not  completely  fill  up  the  impulse  groups  A.  Let  the 
counting  cycle  -  turned  on  by  the  MMV  -  start  just  at 
that  moment  when  the  impulse  group  t  in  counter  2  cor¬ 
responding  exactly  to  a  block  tl  is  over  . 

Case  c)  Let  the  frequency  of  A  be  slightly  larger  than 
in  case  a,  so  that  the  360,000  impulses  of  counter  1 
need  a  bit  more  time  than  is  corresponding  to  an  integral 
number  of  impulses  A.  we  may  have  the  case  that  at 
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counter  2 


counter  1 


counter  2 
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starting  the  counting  according  to  case  a)  we  count  just 
one  additional  impulse  train  t  in  counter  2,  whereas  the 
corresponding  impulse  group  tl  in  counter  1  is  not 
counted  completely. 


Calculation 

To  be  able  to  conduct  calculations  the  impulses  coming 
under  the  period  .  must  be  distributed  equally  to  the 
time  tl.  Therefrom  we  can  calculate  the  lack  (case  b) 
or  the  excess  (casec)  of  impulses  in  counter  2. 

Case  b:  At  first  we  calculate  the  number  of  impulses  n 
that  v?ere  missing  in  counter  2  in  relation  to  counter  1 

n  =  _L  (t  -  t)v 

tl  1  2 

From  the  definition  of  the  impulse  width  t 
t  :  2tj  =  $  :  360 


we  get  with 


A*2  =  n 


the  equation  for  error  2 


Ho  = 


2^  v, 


A<J>2  ...  error  2 

<J>  . phase  position 

Vj  ....  measuring  frequency 

v2  ....  clock  frequency 

a  .  maximum  content  of  counter  1 


For  case  c)  the  sams  calculation  holds,  but  the  error 
is  negative. 
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At  the  device  described  here  the  following  values  were 
fixed  for  a  and  v ^ 

a  .  .  .  .  360,000 
v^  .  .  .  10^  Hz 

From  Eq.  (5.2)  follows  that  this  error  becomes  smallest 
When  either  a  becomes  very  high  or  the  ratio  of  block 
frequency  to  measuring  frequency  decreases. 

Table  5.2  shows  the  values  of  the  two  errors  in  dependence 
on  phase  and  measuring  frequency. 

Table  5.2 


v!  I«z| 

♦  l°l 

M,  H 

M 

103 

3.6 

+0.036 

+0.2 

io3 

180.0 

+0.036 

+0 

io4 

3.6 

+0.36 

+0.02 

104 

180.0 

+0.36 

±° 

io5 

3.6 

+3  •  6 

+0.002 

io5 

180.0 

+3  •  6 

+0 

As  can  be  seen  from  Fig.  5. 

9,  error  2 

is  zero  at  180° 

and  increases 

with  decreasing  phase. 

It  is  to  be  mentioned  that 

the  duration  of  a  counting 

-2 

cycle  is  3.6. 

10  sec;  this 

i  time  can  be  easily  increased 

as  required  by  means  of  the  two  MMV  in  Fig.  5.4. 


CHAPTER  6 


VAT.LEV  MEASUREMENTS 


The  above-ground  registration  of  VLF  signals  showed  the 
maximum  magnetic  field  strength  always  to  run  in  direc¬ 
tion  of  the  mountain  valleys.  This  effect  obviously 
does  not  occur  by  chance  because  it  could  be  proved  in 
all  by-valleys  of  the  Inn-Valley.  A  thorough  explanation 
will  be  given  in  a  particular  thesis. 

In  spring  1972  we  started  preparatory  measurements. 

Several  frequencies  were  considered  (15.1;  50;  128.25; 

278  kHz)  to  be  able  to  recognize  a  possible  frequency 
dependence.  In  fact,  a  new  unexpected  effect  occurred. 

For  the  higher  frequencies  the  maximum  magnetic  field 
strength  was  always  normal  to  the  valley  direction.  Fig.  1 
shows  these  results.  The  contour  line*;  show  the  course  of 
the  valleys  in  the  surroundings  of  Innsbruck  and  the  vec¬ 
tors  show  the  direction  of  the  maximum  magnetic  field 
strength.  If  the  valleys  become  too  narrow  all  vectors 
point  into  valley  direction. 

For  a  hollow  waveguide  the  mode  with  the  greatest  wave¬ 
length  is  a  1 E-mode,  i.e.,  a  mode  with  the  maximum  magne¬ 
tic  field  strengck  following  the  waveguide  axis.  So  we 
tried  to  interprete  the  effects  by  a  modified  waveguide 
theory.  The  only  difficulty  was  to  recognize  the  valleys 
as  hollow  waveguides.  In  several  publications  /23,  24, 

25/,  the  problem  of  waves  over  rough  surfaces  is  solved  by 
expanding  the  irregularities  of  the  surface  in  Fourier- 
series.  The  wave  equations  are  solved  for  sinoidal 
modulated  surfaces  and  then  superponed.  It  is  evident  that 
surface  irregularities  will  induce  new  modes. 


iAftki 


-  62 


4 


In  /24/  we  find  waveguides  of  arbitrary  cross  section 
being  transformed  into  cylindrical  waveguides  by  means 
of  conformal  mapping.  The  problem  now  appeared  rather 
restricted:  The  shape  of  the  valley  has  to  be  transformed 
to  a  cylinder. 


w  *  f(z) 


Fig.  6.2 

The  exact  solution  requires  considerable  numerical 
expenditure.  Besides,  there  is  no  plotter  at  Innsbruck 
University  so  that  the  results  cannot  be  plotted.  For 
the  time  being  we  have  to  be  satisfied  with  an  approxi¬ 
mate  solution  which,  however,  shows  good  physical  clear¬ 
ness. 

The  transformation: 

z-i 

By  means  of  the  function  w  =  — —  the  upper  semi-plane 
is  represented  onto  the  unit-circle  of  the  w-plane 
/18,  26/. 


z-i 

z+i 

Fig.  6.3 


. 


C  3 


We  choose  our  valley  cross  section  In  a  way  that  it  falls 
completely  into  the  upper  semi -plane. 


Fig.  6.4 


Now  the  valley  is  transformed  into  the  interior  of  two 
intersecting  circles  inside  the  unit  circle  which  depend 
on  the  valley  size. 

Waveguides  still  can  be  described  by  the  eigenfunctions 
of  cylindrical  waveguides,  if  the  cross  section  varies 
slightly  from  a  circle  /24/.  In  our  case,  this  means  that 
the  problem  is  solved  principally  when  the  valley  is  re¬ 
presented  onto  a  fairly  circular  area. 


The  representation  function  reads 

z  -  i 

z  +i 


The  circle 

j  z  -  ia| 2  =  a2 


(6.1) 

(6.2) 


in  the  z-plane  is  transformed  into  the  circle 

f  l  } 1  2  4a2 

l  l+2aj  (1  +  2a) 2 


+  v2  = 


(6.3) 


in  the  w-plane. 

For  the  straight  line  Im(z)  =  b  we  obtain  in  the  w-plane 


64 


Considering  e.g.  the  case  a  «  b  «  3,  an  approximately 
circular  representation  of  the  valley  results  in  the 
w-plane .  The  concentric  circles  to  the  circle  K  in  the 
w-plane  are  representing  coordinate  lines  in  the  hollow 
waveguide. 


♦* 


2 


-  C2 


(6.5) 


The  inverse  transformation  of  these  circles  produces 

48  +  49c2 

- - - 

64  -  49c2 

in  the  z-plane.  These  circles  are  shown  in  Pig.  6.5. 


7252c2 


(64  -  49c2)2 


(6.6) 


Fig.  6.5 
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Th«  sane  results  are  obtained  when  using  other  values 
for  a  and  b.  So  we  can  say  that  the  transformation  is 
independent  on  the  size  of  the  valley. 

Conclusion 

It  is  possible  to  describe  a  mountain  valley  by  a 
cylindrical  waveguide-  by  means  of  conformal  mapping. 

From  Fig.  6.5  it  is  to  be  seen  that  near  the  bottom 
the  field  density  is  higher  than  in  the  upper  regions 
of  the  valley. 

An  exact  transformation  of  the  valley  cross  section  into 
a  cylindrical  waveguide  is  possible,  but  requires  a  com¬ 
puter  with  plotter  because  of  the  great  numerical  ex¬ 
penditure. 
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CHAPTER  7 

A  NEW  METHOD  OF  DETERMINING 
THE  ELECTRIC  ROCK  PARAMETERS  ey  AND  o 


7 . 1  General  considerations 

The  procedures  O;  determining  er  and  c  listed  below  are 
based  on  the  dependence  of  the  field  strength  on 

1)  the  dipole  moment  in 

2)  the  distance  R  between  point  of  reception  and 
transmitter 

3)  the  frequency  f 

4)  the  dielectric  constant  er 

5)  the  conductivity  o 

H  »  H(m,  r,  f,  er,  a)  (7.1) 

H  ...  magnetic  field  strength  vector 

To  restrict  this  complexity  it  is  necessary  to  conduct 
the  measurement  in  a  way  that  as  many  of  these  parameters 
as  possible  are  kept  constant. 

Furthermore,  we  have  to  consider  the  problem  whether 
and  how  much  the  individual  parameters  depend  on  each 
other. 

a)  Since  |m|  =  n*I*F  can  always  be  kept  constant  by  read¬ 
justing  the  antenna  current,  the  dependence 

I  in  1  =  |  in  |  (er,  a) 

that  might  possibly  arise  by  attenuating  the  trans¬ 
mitting  antenna  operated  in  resonance,  is  only  of 
theoretical  interest  for  our  case. 


b)  It  is  obvious  that 


er  +  er(m);  a  £  c(m);  £  ^  £(m);  and  R  +  R(m) 

c)  The  question  for  the  dependence 

?  2 

er  =  er(R)  and  o  *  o(R), 

however#  has  to  be  discussed  more  thoroughly  which  is 
to  be  illustrated"B$r  ^teans'“6f"the^  following  figures 
(Figs.  7.1a  and  7.1b). 


Fig.  7.1a 

In  the  first  case  (Fig.  7.1a)  transmitter  and  receiver 
are  so  close  together  that  solely  the  dielectric  con¬ 
stant  Cj  of  the  inhomogeneous  island  in  between  holds 
responsible  for  the  field  at  the  point  of  reception.  If 
the  measuring  distance  is  enlarged  the  surrounding  medium 
characterized  by  e2  becomes  more  and  more  effective. 

The  resulting  dielectric  constant  eres  influencing  the 
field  of  reception  is  composed  of  e l  and  e2.  Fig.  7.1b 
will  illustrate  this  situation.  It  was  on  purpose  that 
all  disturbing  bodies  were  drawn  with  equal  magnitude  and 
shape  to  make  clear  then  even  then  eres  is  depending  on 
the  distance.  In  reality  these  islands  will  be  of  different 
size  and  have  a  different  dielectric^ constant,  too. 
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The  distance  dependence  <cr  er  and  a  (of  course,  we  can 
reflec*:  in  the  same  way  on  the  conductivity)  can  thus 
be  regarded  as  a  fact.  The  dielectric  constant  (or  con¬ 
ductivity,  respectively)  determined  over  a  certain 
measuring  distance  is  a  mean  value  of  all  dielectric  con¬ 
stants  (or  conductivities)  of  this  region. 

Vt  is  therefore  senseless  to  compare  specific  values  of 
dielectric  constant  (or  conductivity) ,  unless  they  are 
related  to  a  certain  distance.  Consequently,  listed 
values  given  by  different  authors  can  only  be  regarded 
as  approximations  and  it  is  necessary  to  compare  measuring 
methods  and  conditions,  to 


d)  The  dependence 

?  ? 
e  =  tr(r)  and  a  -  o(f) 

is  of  ••rest  importance  witn  regard  to  the  two-frequency 
method  discussed  below.  The  frequency  dependence  of  er 
and  the  consequent  dependence  of  the  index  of  refraction 
on  the  frequency  is  called  "dispersion”  in  optics.  In 
this  case,  however,  wavelengths  involved  are  so  small  th.t 
they  can  be  compared  wl th  molecules  and  resonance  vibra¬ 
tions  occur.  At  low  frequencies  there  are  no  such  vibra¬ 
tions;  the  approximate  limit  lies  at  about  30  MHz 
U0  =  10  m,  SW-band) ,  Eelow  this  frequency, calculations 
are  performed  with  the  static  value  of  the  dielectric 
constant. 

In  contrast  to  this  assumption  a  certain  frequency  re¬ 
sponse  was  observed  in  the  lower  frequency  range 
(f  <  100  kHz)  owing  to  the  so-called  partition  polariza¬ 
tion.  Among  others,  STRUTT  /27/  ard  recently  WtJBKING  /28/ 
and  GRISSENANN  / 2 9/  have  observed  this  phenomenon.  In 
/ 2 8 /  and  /29/,  the  (cf,f)  and  the  (a ,f) -diagrams  show  a 
flattening  of  the  curves  towards  higher  frequencies  (f >  100  kHz) 
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Whether  or  not  the  two-frequency  method  is  admissible 
for  the  determination  of  er  and  a  depends  on  the  above 
conclusions. 

For  a  comparison  let  us  briefly  refer  to  methods  that 
have  already  been  worked  out  by  co-workers  of  the  VLF- 
Project: 


7.2  Method  according  to  NESSLER  /10/ 

If  the  relative  value  G  (to  be  defined  below)  is  plotted 
as  ordinate  for  a  definite  frequency  and  a  definite 
measuring  distance  every  point  in  the  diagram  will  con¬ 
stitute  a  pair  of  values  (er,  o)  belonging  together.  To 
obtain  H90  we  must  know  exactly  the  properties  of  the  re¬ 
ceiving  equipment , and  the  calibration  factor  establishing 
the  interrelation  between  actual  field  strength  and 
measured  voltage  must  be  determined  as  accurately  as 
possible.  Since  the  calibration  factor  comprises  quite 
a  number  of  measuring  values,  the  uncertainty  factor  is 
rather  high.  The  resolving  power  of  the  method  is,  how¬ 
ever,  good  even  at  small  measuring  distances  of  about  20  m. 


7.3  Method  according  to  KELLNER  / 1 4 / 

Based  on  an  adequate  representation  of  the  formulas  for 
the  vector  components  of  the  field  strength,  KELLNER  /14/ 
found  a  possibility  to  determine  the  complex  wave  number 
from  relative  values  and  thereby  calculate  er  and  a. 

The  main  disadvantage  of  this  method  is  its  failure  in 
the  inhomogeneous  measuring  site.  According  to  the  author, 
the  method  can  only  be  used  up  to  a  frequency  of  about 
350  kHz  and  a  measuring  distance  of  50  m.  This  is  due  to 
the  fact  that  with  decreasing  wave  length  the  inhomogeni- 
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ties  form  an  ever  increasing  hindrance  for  the  field  j 

» 

which  thus  is  more  and  more  distorted.  It  is  chiefly  the  1 

bearing  angle  ♦  indicating  the  direction  of  maximum  field  \ 

strength  and  playing  a  dominating  role  in  the  evaluation  j 

l 

formulas  that  is  distorted,  thus  heavily  falsifying  the  j 

result.  i 


7.4  New  Method  according  to  KLEIN  /17/ 

Two-frequency  method 

To  be  able  to  make  quantitative  statements  on  the  beha¬ 
vior  of  electromagnetic  waves  penetrating  zones  of  in- 
homogenity  it  is  necessary  to  know  the  electric  rock 
parameters  er  and  o.  Real  and  imaginary  parts  of  the  com¬ 
plex  wave  number  k*  and  of  the  complex  refraction  index 
n*  can  be  calculated  from  these  parameters.  Therefore  we 
had  to  find  a  method  delivering  workable  results  even  in 
an  inhomogeneous  medium. 

The  following  requirements  had  to  be  fulfilled: 

1)  The  method  must  be  practicable  in  inhomogeneous  medium. 

2)  In  its  formulas  and  diagrams  it  should  contain  only 
quotients  of  absolute  measuring  values. 

A  favorable  application  of  a  possible  method  made  us 
request  two  more  conditions: 

3)  Measurements  should  be  able  to  be  conducted  easily 
and  quickly. 

4)  The  values  er  and  o  to  be  determined  should  possibly 
be  readable  at  the  measuring  site  thus  enabling  a 
short-time  change  in  the  measuring  program. 

A  solution  of  these  problems  was  found  using  the  field 
strength  ratio  G  (as  was  done  with  the  two  methods 
described  above)  considering,  however,  two  different 
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measuring  frequencies.  According  to  /10/  G  is  defined 
as 

Hn  Hp 

G  =  r-fl-  *  rr“  tan  0  (7.2) 

m90  "0 

HQ  ...  field  strength  maximum  at  6  =  0° 

H90  . .  field  strength  maximum  at  e  =  90° 

HR  . . .  field  strength  in  the  direction  towards  the 

transmitting  antenna 
Hq  ...  field  strength  vertical  to  HR 
O  ...  angle  between  transmitting  antenna  and 

radius  vector 

-¥ 

Since  the  field  strength  vector  H  and  thus  all  its  com¬ 
ponents  depend  on  frequency,  dielectric  constant,  con¬ 
ductivity,  and  distance  to  the  transmitter 

G  e  G(f ,  ^0 

Owincr  to  the  above  mentioned  reasons  the  measuring  distance 
transmitter  -  point  of  reception  must  be  kept  constant,  i.e. 

G  =  G(f ,  c ^ ,  a) 

The  above  mentioned  frequency  independence  of  er  and  o 
at  one  and  the  same  measuring  distance  reduces  the  de¬ 
pendence  of  G  to 

G  =  G(f ) 

at  er  =  const,  and  o  =  const. 

If  groups  of  curves  with  the  parameters  er  and  o,  resp.  , 
are  plotted  in  a  diagram  with  the  ordinate  G  =  G(er,o) 
at  a  frequency  fj  and  with  the  abscissa  G  =  G(er,a)  at 
a  frequency  f2 ,  every  point  of  intersection  of  the 
er-parameter  group  with  the  a -parameter  group  constitutes 
a  (er,o)-pair  of  values  belonging  definitely  together. 
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For  our  measurements#  the  following  values  were  chosen: 
fj  =  500  kHz;  f2  ■  1000  kHz;  R  =100  m  (50  m),  thereby 
taking  into  consideration  the  geometry  of  the  gallery# 
the  resolving  power  at  a  definite  measuring  distance 
(splitting  of  the  curves  of  the  diagram) #  and  the  power 
of  the  transmitter.  For  G(500)  (i.e.  G  at  500  kHz)  and 
G(1000)  (i.e.  G  at  1000  kHz)  determined  by  measurement# 
dielectric  constant  and  conductivity  can  be  taken  from 
Fig.  7.2.  For  other  measuring  distances  and  frequencies 
corresponding  diagrams  have  to  be  plotted. 


antenna  is 
the  determi- 
values  of  e. 


thus  considerably  reducing  measuring  expenditure# 
because  only  the  frequency  has  to  be  interchanged. 

c)  Using  the  relation 


the  transmitting  antenna  is  first  adjusted  to  6  =  0° 
and  then  to  e  »  90°.  At  both  positions  the  maximum 
field  strength  is  measured  for  both  frequencies. 

Out  of  these  three  possibilities  it  is  the  third  one  that 
offers  the  decisive  advantage  in  comparison  to  the  methods 
by  NESSLER  and  KELLNER. 


There  are  three  possibilities  to  measure  G: 
a)  With  the  relation 

G  =  —  tan  e 

He 

any  angle  position  #  of  the  transmitting 
admissible.  For  e  90°,  tane  -►  «,  i.e., 
nation  of  G  becomes  inaccurate  for  large 

bf  If  »  45° 

hR 

G  =  - 

He 
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Both  for  6=0°  and  for  6  =  90°  the  field  is  linearly 
polarized.  When  measuring,  the  H0  and  H90  directions  are 
sought  and  the  amount  is  measured.  That  is  why  a  devia¬ 
tion  owing  to  disturbances  can  by  no  means  influence  the 
accuracy  of  the  value  G  -  Hq/H50. 


7.5  Practical  application  of  the  new  method 

Numerous  measurements  were  executed  in  the  mine  of  St.Ger- 
traudi.  Whereas  for  most  measuring  sites  a  homogeneous 
rock  structure  could  be  assumed,  the  measurements  in  "Mor- 
genschlag"  produced  (cf.  Fig.  7.3)  surprisingly  different 
values;  results  are  listed  in  Table  1. 


Table  Is  R  «  100  m 


exp. 

no. 

( 

measuring 

distance 

f 

[kHz] 

[jvJ 

U90 

[»v] 

G(500) 

GC1000) 

er 

0 

[S/m] 

97 

M8.2  -Ml 2 

500 

1000 

290 

78 

880 

420 

0.33 

0.19 

20 

9  •  10”4 

101 

Ml 2  -M8.2 

500 

1000 

505 

165 

1300 

810 

0.39 

0.21 

20 

5  •  10~4 

100 

Mil.  3  -M9 

500 

1000 

490 

108 

790 

390 

0.62 

0.28 

12 

io“4 

104 

M9  -  Mil. 3 

500 

1000 

485 

120 

780 

400 

0.62 

0.30 

10 

10“4 

99 

Mil.  3-MU 

500 

1000 

125 

55 

138 

143 

0.91 

0.39 

8 

6 • 10~5 

103 

MU-MU. 3 

500 

1000 

445 

120 

610 

360 

0.73 

0.33 

8 

8  •  10*”5 

98 

MU.3-M12 

500 

1000 

530 

92 

510 

205 

1.04 

0.45 

4 

2 • 10-5 

102 

i  M12-MU.3 

! 

500 

1000 

480 

85 

510 

210 

0.94  j 
0.41 

6 

5 • 10-5 
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The  number  of  experiment  refers  to  Fig.  7.3.  On  com¬ 
paring  exp.  no.  97  with  exp*  98  perpendicular  thereto 
we  realize  that  the  two  dielectric  constants  are  5  :  1 , 
the  conductivities ,  however,  nearly  50  s  1.  From  these 
remarkable  differences  we  can  draw  the  conclusion  that 
a  disturbance  runs  through  the  medium.  The  recognition 
of  such  disturbances  means  a  real  advantage  for  geologi¬ 
cal  prospecting.  The  method  described  above  together  with 
the  method  of  images  thoroughly  discussed  in  /17/  has 
already  produced  promising  results. 


7.6  A  further  possibility  to  locate  disturbing 
zones  in  rock 


The  following  characteristic  resulted  from  the  fore- 
mentioned  measurements  for  the  determination  of  the 
electric  rock  parameters  er  and  a  in  "Morgenschlag" : 

At  the  execution  of  measurement  106  (78)  (cf.  Fig.  7.3) 
from  M8.2  to  M12  the  direction  of  maximum  field  strength 
(major  axis  of  receiving  ellipse)  rotates  by  180°  together 
with  the  rotation  of  the  transmitting  antenna  whereas  at 
measurements  perpendicular  thereto  the  maximum  field 
strength  shows  some  kind  of  pendulous  motion.  The  rota¬ 
tion  occurring  in  one  direction  up  to  a  certain  angle  of 
rotation  e  of  the  transmitting  antenna  suddenly  reverses 
so  that  the  receiving  antenna  tilts  back  to  the  starting 
position  without  having  done  a  l80°-rotation. 

Fig.  7.4  illustrates  quite  clearly  this  extraordinary 
behavior.  The  curves  do  not  proceed  -  as  could  be  expected 
from  theory  -  from  the  first  to  the  third  quadrant,  but 
from  the  first  to  the  second  one  at  the  forward  measure¬ 
ment  and  from  the  fourth  to  the  third  quadrant  at  the  back 
ward  measurement  (with  exchanged  transmitter  and  receiver) 
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The  curve  of  the  backward  measurement  is,  as  it  were, 
the  "image"  of  the  forward  measurement  curve. 

An  exact,  purely  mathematical  description  of  this  effect 
would  be  extremely  difficult,  because  too  many  unknown 
quantities  or  variables  are  involved.  Certainly,  numerous 
papers  were  published  on  this  subject  but  finally  calcula 
tions  are  always  based  on  the  assumption  of  quite  special 
facts.  The  main  object  of  our  method,  however,  was  to 
give  at  least  a  qualitative  explanation  for  the  "pendu¬ 
lous  motion"  by  means  of  a  clearly  arranged  geometrical 
construction  using  the  method  of  images. 

The  development  of  the  theory  was  stimulated  by  /37, 

38,  39,  40/. 


7.7  Theory  of  images 

If  a  quiescent  charge  is  situated  in  direct  neighborhood 
of  a  conducting  infinite  surface  a  field  distribution 
according  to  Fig.  7.5a  results.  This  consideration  shall 
app  f  analogously  for  moved  charges  (Fig.  7.5b).  The 
re1  cion  between  the  magnetic  dipole  moment  and  its 
.  age  is  shown  in  Fig.  7.5c. 

There  are  two  measuring  possibilities  in  the  neighbor¬ 
hood  of  a  disturbing  zone: 

Type  1 :  Measurement  approximately  parallel  to  disturbing 


zone  , 

Jn4o' 


Fig.  7.6 
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Let  SA'  be  the  mirror  image  of  SA.  At  EA,  two  signals 
are  received  superponing  one  another.  The  two  field 
vectors  H  and  H'  generally  having  different  directions 
form  a  resulting  vector  Hres  with  another  direction. 

Type  2  s  The  conducting  zone  is  situated  between  trans- 
mitting  antenna  and  receiver 


Here,  the  image  SA'  develops  from  the  charge  distribu¬ 
tion  induced  by  SA  in  the  conducting  layer  and  acting  as 
an  image  of  SA  at  SA'.  The  field  strength  at  EA  is  again 
composed  of  the  signal  originating  at  SA  and  SA'. 

The  good  agreement  between  theory  and  experiment  shall  be 
shown  by  means  of  two  examples: 

Pig.  7.8  is  an  example  for  type  1.  The  measuring  distance 
is  parallel  to  the  disturbance  determined  in  "Morgenschlag" . 

Pig.  7.9  represents  type  2  with  the  disturbance  between 
transmitter  and  receiver. 

Prom  the^  above  figures  we  recognize  that  the  auxiliary 
vector  H'  explains  rather  well  the  deviating  direction 
of  Hex. 

If,  however,  H  and  H'  include  an  angle  of  approximately 


FIG.  7.5 


180°  some  difficulties  in  interpretation  will  arise. 

-4 

The  direction  of  the  resulting  vector  Hex  coincides  with 

-4  4  -4 

H  or  H'  at  type  1  (exp.  no.  106),  whereas  Hex  is  perpendi 

-4  -4 

cular  to  H  and  H'  at  type  2  (exp.  no.  88) .  So  we  have  to 

-4  4 

assume  that  a  phase  reversal  in  respect  to  II 1  (or  H) 
occurs  at  the  second  case,  thereby  giving  a  good  explana¬ 
tion  for  the  pendulous  motion  characterized  by  the  re¬ 
versal*  of  rotation . 

The  method  developed  in  the  mine  of  St.  Gertraudi  was 
already  successfully  applied  in  other  mines  (Schwaz/Tyrol 
Bleiberg/Carinthia)  thereby  confirming  our  theory.  A  com¬ 
mercial  application  of  this  method  for  ore  prospecting, 
however,  would  require  further  experiments  to  investigate 
the  influence  of  watercourses,  rails,  and  other  technical 
equipment . 
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CHAPTER  8 

ENTERING  MECHANISM  OF  ELECTROMAGNETIC  WAVES 
INTO  A  DISSIPATIVE  MEDIUM 


8 . I  Introduction 

In  December  1972  the  Investigations  on  this  subject  were 
finished.  The  basic  theory  which  will  be  briefly  referred 
to  was  already  published  in  Final  Report  1971  together 
with  the  development  of  an  autorange  receiver. 

Problem; 

Measurements  /ll/  of  the  broadcasting  transmitter  O-Re- 
gionax  along  a  gallery  of  the  mine  St.  Gertraudi  showed 
a  local  period  of  the  maximum  magnetic  field  strength 
(major  axis  of  a  rotary  field  ellipse) .  The  transmitter 
(629  kHz)  is  situated  at  a  distance  of  30km  from  the  mine. 

Repeated  measurements  confirmed  this  unexpected  result 
which  was  tried  to  be  explained  / 30/  by  developing  a 
theory  for  an  electromagnetic  wave  with  oblique  inci¬ 
dence  on  a  boundary  surface. 


8.2  Equations  of  the  magnetic  field  strength  in  medium 

The  equations  used  by  Schumann  /30/  are  specified  in  /ll/. 
For  H  polarized  in  the  direction  of  the  plane  of  inci¬ 
dence  (in  air)  they  take  the  form 

Ez(x,y)  =  Ege  exp{ jut  +  f(x,y)} 

cos$ 


H  (xfy)  =  -E®  — L  exp{ jut  +  f (x,y) } 


Hy  (xfy) 


wi 

*  -Ese  expCjut  +  f  (x,y) } 


(8.1) 


Wi 
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for  H  polarized  vertically  to  the  plane  of  incidence 
(in  air)  they  read: 

Hz(x,y)  *  Hge  exptjwt  +  f(x,y)} 

Ex(x,y)  »  Hse*W1  cos^j  expCjwt  +  f(x,y))  (8.2) 

Ey(x,y)  =»  Hse*wi  Bln#!  exp{j*»t  +  f(x,y)} 

Ese  . . .  electric  field  strength  of  the  incoming  wave 
Hg«  . . .  magnetic  field  strength  of  the  incoming  wave 
W1  . ...  v *ave  resistance  in  air  ■  377  0 
* l  ....  angle  of  incidence 
. . circular  frequency 


With  the  boundary  conditions  the  complex  law  of  Snellius 
is  obtained 

Wi 


sin* 


l 


W, 


sin*. 


W2  =  pw  exp 


o/2 


Pw, 


w/y2e2  4 

- - /I  + 


(8.3) 


c0  Aj2  +  w2e  2  c2 


“C0C2; 


arctan- 


U)C  0^2 


w2  ....  wave  resistance  in  medium  (rock) 
£ 2  ....  dielectric  constant  in  medium 
e0  ....  absolute  dielectric  constant 
u2  ....  permeability  in  medium  (rock) 
a  ....  conductivity  in  medium  (rock) 
cQ  ....  vacuum  light  velocity 


•  nn  iiaii  iinr  uriTi  Ban  eSi  nnfi  r  nt  ' '  r 
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From  the  boundary  conditions  we  can  conclude  the  com¬ 
plex  transmission  coefficients  for  the  case  of  parallel 
polarization  (1+Tp)  and  vertical  polarization  (1  +  rv) , 
respectively 


After  laborious  calculations  the  following  equations  for 
the  magnetic  field  strength  are  obtained 

Hxd  =  pi+r  >'p  *W1  *Hp*e^*exp{  j  |wt  +  ?f  +  (Ip  -  e/2  |  } 

P 

Hyd  “  Pi+rp*8in,?i  ‘Hp-e^expC  j  |ut  +tt  +  Op  1  >  (8.6) 

Hzd  "  Pi+r^H^-e^exptjUt  +  Cp  -  cp  +  cv  -  a | } 


p,,r  ...  amount  of  the  transmission  coefficient  of 

1  *  p  -> 

the  electric  field  strength  with  H  polarized 
parallel  to  the  plane  of  incidence 


pi+rv 


**♦ (x#y) 

+  -> 

0p=-k x . . 

k . 

Hg/Hg0  .  .  .. 

tp >Cg  • • • 

n . 


amount  of  the  transmission  coefficient  of 

-V 

the  magnetic  field  strength  with  R  polarized 

vertically  to  the  plane  of  incidence 

attenuation  factor 

position  function  of  the  phase 

vector  of  wave  number 

amplitudes  of  the  magnetic  field  strength 
of  the  incoming  elliptically  polarized  wave 
phases  of  the  complex  transmission  coef¬ 
ficients 

phase  shift  of  the  incoming  elliptically 
polarized  field 


Furthermore , 


(A) 


5  r^J-  a 


is  valid  with  high  accuracy  (error  less  than  1  %) . 


In  the  following  paragraph  the  reasons  leading  to  a 
local  periodicity  are  dealt  with. 


8.3  Reasons  for  the  local  periodicity 
Assumption; 

The  rock  is  no  longer  assumed  to  be  homogeneous  but 
bearing  a  reflecting  layer  which  is  shown  in  Fig.  3a. 

By  superposition  of  the  incoming  and  the  reflected  wave 
a  local  periodicity  is  generated.  Further  calculations 
show  the  possibility  of  determining  the  angle  between  re¬ 
flected  layer  and  gallery  axis  from  the  three  occurring 
modes; 

a)  k^.  in  the  first  quadrant  of  the  complex 

period  only  in  z ' -direction.  Since  kxi>r  has  the  same 
direction  as  kxi  no  periodicity  can  occur  in  x‘ -direction. 
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b)  kr  in  the  second  quadrant 

x1-  and  z '  -periodicity .  and  k  «re  anti- 

parallel  to  k*!  and  k2..  ' 

c)  kr  in  the  third  quadrant 

period  only  in  x* -direction.  k*.^  antiparallel  to 
^St'»  ^*'rr  Parallel  to  kji. 

(Th.  dashed  quantities  relate  to  the  <S5idinite  system 
determined  by  the  gallery  system.) 

Fig.  3b  shows  the  possible  modes. 

The  following  relations  between  the  angles  can  be  read 
from  Pig.  3a: 


a)  0<K*i+2Y-e<i 

b)  <><**£«  e-2y<- 

2  2 

C)  0  <  K  «=  2y  -  (e  +  i)  <  i 

2  2 


(8.7) 


8.4  Determination  of  v 


From 


8p  -  -kx 


kx»  and  k2i  can  be  obtained. 


_1_  J  cos24i 


(8.8) 


V 


377  —  .  S 

cn  *2  Wtnt) 


(8.9) 

T  ”**  sloPe  angle  o  . rock  conductivity 

angle  of  incidence  c  ...  dielectric  constant  of 

the  rock 


Medium  1  (Luff) 


leren 


Wellentypen.die  bei  Anwesenheit  einer 
reflekt/erenden  Schicht  im  Medium  2 
dort  auftreten  konnen 
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From 

k,.,  =  |  k  |  costc 

kz*,r  *  |  k|  sinie 

and 

|k|  =  /  kx  i2  +  k2i2 

we  can  calculate  both  the  period  in  x' -direction 


(8.10) 


Px’ 


-2* 


kx*  +  kx' ,r 
and  the  period  in  z' -direct ion 


(8.11) 


(3.12) 


Thus,  it  is  possible  to  determine  kx,^r  and  kz»^r  from 
the  rock  parameters,  and  the  period  k  can  be  determined  by 
means  of  Eq.  (8.10)  and  finally  y  is  calculated  from 
Eq.  (8.7)  if  the  mode  is  known. 

Measurements  showed  a  periodicity  in  z'  thus  proving 
the  occurrence  of  case  b) . 

For  e  *  10,  o  =  10" 3  S/mm,  f  =  629  kHz  and  ^  84° 


y  =  36,7° 


The  above  described  method  fails  if  wavelength  and  dimen¬ 
sion  of  the  reflecting  surface  are  in  the  same  order  of 
magnitude.  This  fact  could  be  proved  by  experiments  at  low 
frequencies  (278  kHz).  Eqs.  (8.8)  and  (8.9)  show  that 
kx»  and  kz»  are  proportional  to  u,  the  periods  Px*  and 
Pz»  being  inversely  proportional  to  u. 

Fig.  4b  is  a  schematical  representation  of  the  surroundings 
of  the  measuring  site  in  the  mine  of  Grosskogel. 
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Theory  says  that  for  plane  waves  the  x-  and  u -waves  are 
not  attenuated,  whereas  the  y-  and  v-waves  are.  In¬ 
fluencing  the  test  results  by  the  v-wave  is  impossible 
on  account  of  the  large  distance. 


8.5  Transmission  coefficients 


The  amounts  of  the  transmission  coefficients  were  com¬ 
puted  for  different  values  of  a  and  represented  in  de¬ 
pendence  on  ♦  ;  in  the  diagrams  4c  and  4d. 


(p1+r^>el  always  remains  <1,  whereas  (o1+r  can 

grow  much  larger  them  1. 


The  present  calculations,  however,  are  not  contradictory 
to  the  law  of  energy  preservation,  since  the  transmission 
coefficient  of  the  pointing  vector  S  is  smaller  than  1 
in  any  case  as  is  shown  in  Fig.  4e. 


A  possible  explanation  can  be:  The  electric  field  in  the 
dissipative  medium  is  immediately  transformed  to  an 
eddy  current  field.  Thus,  a  transformation  of  electric 
into  magnetic  energy  occurs.  This  is  in  good  agreement 
with  the  unexpectedly  high  measuring  values  in  winter. 


On  the  one  hand,  the  boundary  surfaces  are  covered  by  a 
well  conducting  snow  layer,  on  the  other  hand,  the 
attenuation  factor  diminishes  on  account  of  the  de¬ 
creased  humidity  in  winter. 


8.6  The  pointing  vector  S 

S  was  investigated  for  H  polarized  parallel  to  the  plane 
of  incidence  and  was  found  to  be  by  no  means  stationary 
any  more.  The  top  of  the  vector  S  is  following  an  ellipse 
beiner  situated  in  the  plane  of  incidence. 


<g.  4d 


Richtung  ,die 


zum  Einf allslot  efnnimmt 


Fig.  4g 
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For  the  general  case  of  an  elliptically  polarized  wave 
S  is  describing  an  oblique  elliptical  cone  in  the 
dissipative  medium.  S  is  complex;  the  mean  time  values 
of  ReSfc  and  ImS t  were  calculated 


So  only  half  of  the  total  energy  is  measurable. 

Finally ,  the  direction  of  ReSfc  in  respect  to  the  vertical 
incidence  was  examined  and  represen ted^ir  Fig .  4g  for 
different  a  in  dependence  on  +. . 

I 

I 
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CHAPTER  9 


ABOVE-GROUND  AND  SUBSURFACE  LONG-TIME 
REGISTRATION  OF  VLF  TRANSMITTERS 


Ba’&ed  4on  some  preparatory  experiments  of  JH.»_  GRADL  per¬ 
formed  in  1965,  the  wellknovm  British  VLF  transmitters 
GBR  at  16.0  kHz  and  GBZ  at  19.6  kHz  were  registered. 

The  main  purpose  of  this  registration  is  to  determine 
whether  or  not  a  surface  wave  of  measurable  amplitude 
is  existing.  Subsurface  measurements  shall  make  it 
possible  to  eliminate  above-ground  atmospherics. 

To  measure  the  GBR  and  GBZ  signals  a  tube-operated  se¬ 
lective  levelmeter  with  a  total  band  width  of  80  Hz, 
which  is  just  including  the  frequency  jump  of  50  Hz  at 
frequency  modulation  was  used.  The  simple  rectification 
used  at  the  beginning  proved  to  be  insufficient  because 
the  noise  amplitudes  amounted  to  a  multiple  of  the  ef¬ 
fective  amplitudes  at  times  of  greater  atmospherics 
(thunderstorms) .  Only  a  special  decoder  in  two  stages 
with  time  constants  of  0.1  and  10  sec  allowed  an  un¬ 
disturbed  reception  of  signals. 

The  long  integration  time  constant  has  the  disadvantage 
of  delivering  a  distorted  pattern  if  the  signal  is 
amplitude-keyed,  the  distortion  being  due  to  the  time 
intervals  (in  the  order  of  seconds) .  There  are  always 
2  hours  of  information  followed  by  1  hour  of  outage 
which  showed  up  to  be  a  negligible  loss  of  information. 

Since  the  tube-operated  receiver  did  not  maintain  the 
adjusted  receiving  frequency  with  sufficient  accuracy 
for  a  period  of  several  days  (Owinn  to  unevitable  tempe 
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rature  fluctuations)  a  varicap  was  added  to  the  oscilla¬ 
tor  thus  providing  the  possibility  of  wobbling  a  certain 
frequency  band  externally.  Fig.  9.1  shows  the  block 
diagram  of  the  complete  equipment. 

Long  ferrite  rod  antennas  which  were  optimized  ade¬ 
quately  have  proved  successful  /9,  31,  32/.  The  essential 
piece  of  the  equipment  is  the  selective  levelmeter  with 
an  input  sensitivity  of  -ll  neper  corresponding  to 
about  iO  mV.  The  saw-tooth  voltage  of  a  Wavetek  142 
with  10  min  time  constant  was  used  as  sweep  control. 

The  output  signal  of  the  levelmeter  was  amplified  by 
14  dB  to  max.  4  V  by  a  PM  5175  (Philips)  to  keep  down 
the  error  influence  of  the  threshold  voltage  of  the 
decoding  diodes.  The  decoded  signal  is  registered  on  a 
line  recorder  and  constitutes  a  spectral  distribution 
of  the  VLF  signals  in  the  requested  frequency  range. 

Fig.  9.2  shows  a  par  of  such  a  registration  .ape. 

One  measuring  point  per  transmitter  for  every  10  minutes 
proved  just  sufficient,  an  increase  in  wobbling  fre¬ 
quency  would  have  caused  a  decrease  of  the  registered 
maximum  amplitudes  on  account  of  the  integration  time 
constant. 

The  signals  of  all  receivable  VLF  transmit  ars  in  the 
world  within  the  frequency  range  of  14.5  to  20.5  kHz 
were  received  above  and  below  the  earth  surface.  The 
object  was  to  find  out  whether  over  a  propagation  di¬ 
stance  of  1300  km,  several  100  m  below  the  surface,  the 
sky  waves,  i.e.  the  signal,  that  is  reflected  by  the 
ionosphere,  form  the  main  part  of  the  received  signal 
or  if  the  surface  waves  represent  a  considerable  amount 
of  signal  or  the  signal  at  all. 

For  this  purpose  the  course  of  the  signal  amplitude  of 
the  GBR  transmitter  (16,0  kHz)  was  plotted  in  the  same 
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scale  during  24  hours  from  both  a  measurement  below 
and  above  the  earth  surface  (Fig.  9.3).  The  two  curves 
were  measured  on  two  different  days,  because  only  one 
reception  apparatus  was  available.  A  determination  of 
a  mean  value  of  amplitudes  over  a  period  of  several 
days  was  impossible ,  since  this  would  have  blurred  the 
characteristic  sudden  increase  of  ampli ti  h-  at  sunrise. 

The  two  curves  may  be  regarded  as  standing  tor  a  n-.-'-l  w 
oi  curves  plotted  under  the  pars  condition;;,  From  * he 
curves  it  is  obvious  that  the  wave  in  soil  t curve  2) 
is  certainly  the  sky  wave,  because  the  attenuation 
measured  at  the  very  constant  amplitude  during  daytime 
corresponds  exactly  to  the  precalculated  attenuation 
of  the  sky  wave.  Furthermore,  the  course  of  amplitude 
at  sunrise  and  sunset  characteristic  for  the  sky  wave 
can  be  determined  equally  above  and  below  the  surface. 

If  the  fluctuating  night  course  of  several  measuring 
series  below  and  above  the  surface  is  averaged,  the 
same  amplitude  ratio  as  during  daytime  can  be  calcu¬ 
lated.  If  there  were  a  signal  portion  of  sinnificant 
magnitude  that  is  independent  on  the  ionosohere,  the 
amplitude  ratio  would  char  :,v  r,:\  .  count  of  the  different 
attenuation  for  surface  and  sky  waves.  Thus,  it  is  evident 
that  the  ionosphere  has  a  significant  influence  on  the 
diurnal  course  of  amplitude  even  at  subsurface  reception. 

The  dashed  parts  of  the  curve  in  Fia.  9.3  are  estimated 
values  corresponding  to  the  AM-periods  of  the  transmitter. 
This  characteristic  course  could  be  observed  during  all 
the  time  of  registration  (several  months)  in  practically 
the  same  way.  It  must  be  mentioned  that  at  night  the 
amplitude  is  considerably  smaller  than  at  daytime. 

Plotting  the  GBZ  transmitter  at  19.6  kHz  with  the  same 
location  brought  about  quite  different  results  (Fig.  9.4). 


Sagezahn -  _ ry  akustische 

generator  U L,  Kontrotle 
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In  Fig.  9.4  the  night  course  between  7  p.m.  and  7  a.m. 
is  plotted.  The  diurnal  course  of  the  two  transmitters 
is  approximately  identical  (GBZ  amplitude  is  by  about 
20  %  smaller  than  GBR  amplitude)  whereas  the  GF?  night 
amplitude  comes  up  to  the  value  of  the  day  amplitude. 

In  some  nights  the  day  values  were  even  exceeded. 

The  upper  diagram  of  Fig.  9.4  shows  the  characteristic 
conditions  of  the  above-ground  registration,  the  diagram 
below  shows  the  interpretation  of  a  subsurface  registra¬ 
tion. 

Since  there  was  only  one  device  available  to  plot  all 
these  diagrams  an  exact  coincidence  of  the  reception 
characteristics  is  impossible.  Two  comparable  diagrams 
had  to  be  measured  on  two  different  days. 

To  eliminate  accidental  fluctuations,  a  greater  number 
of  above-ground  and  subsurface  diagrams  was  interpreted 
always  showing  this  distinct  difference  between  the  two 
neighboring  frequencies  of  16  and  19.6  kHz.  It  is 
principally  possible  to  take  an  average  of  the  different 
similar  diagrams  and  would  equalize  the  slight  fluctua¬ 
tions  of  the  night  course,  but  the  characteristic  sunrise 
peaks  would  also  be  blurred. 

For  a  further  examination  of  these  effects  it  is  abso¬ 
lutely  necessary  to  register  the  signal  amplitudes 
above  and  below  the  earth  surface  simultaneously.  Two 
battery-operated  registration  devices  are  already  under 
construction,  using  a  magnetic  tape  instead  of  the  line 
recorder  as  a  data  carrier.  Thereby  it  will  be  possible 
to  show  up  clearly  one  of  the  expected  advantages  of  a 
subsurface  registration,  namely  the  fact,  that  artificial 
or  man-made  noise  peaks  are  more  heavily  attenuated  than 
the  measuring  signal  penetrating  from  outside.  This 
effect  was  partly  determined  with  the  simple  equipment 


90  - 


already  and  is  due  to  the  fact  that  the  signal  of  the 
tw«  transmitters  arrives  as  a  farfield  at  the  point  of 
reception  whereas  the  especially  disturbing  noise  peaks 
from  the  near  neighborhood  arrive  as  a  near field  with 
a  multiple  of  the  measuring  amplitude.  This  is  explained 
by  the  following  two  diagrams  (Figs.  9.5  and  9.6). 

Fig.  9.5  shows  an  above-ground  registration  during  a 
local  thunderstorm;  it  can  be  seen  that  instead  of  the 
zero  line  there  is  a  beam  of  numerous  noise  peaks  which 
come  through  very  distinctly  in  spite  of  the  efficient 
noise  inverter. 

Fig.  9.6  represents  the  equivalent  subsurface  registra¬ 
tion  showing  an  absolutely  undisturbed  diagram  when 
pa  sing  frequency  ranges  without  transmitting  signals. 

An  exact  correlation  of  these  diagrams  is  again  im¬ 
possible,  as  they  were  plotted  at  different  times. 
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